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We estimate the signal-to-noise ratios (SNRs) that one would expect to measure from coalescing bi¬ 
nary black hole (BBH) systems for the following three broadband gravitational-wave observatories: 
initial and advanced ground-based interferometers (LIGO/VIRGO) and space-based interferome¬ 
ters (LISA). We focus particularly on the highly relativistic and nonlinear merger portion of the 
gravitational-wave signal, which comes after the adiabatic inspiral portion and before the ringdown 
portion due to the quasinormal ringing of the Hnal Kerr black hole. 

Ground-based interferometers can do moderate SNR (a few tens), moderate accuracy studies 
of the dynamics of merging black holes in the mass range (a few)M 0 to ~ 2 OOOM 0 . LISA, by 
contrast, can do high SNR (a few x 10^), high-accuracy studies of BBH systems in the mass range 
10® M 0 < {1 z)M < 10® M 0 , where 2 is the binaries’ cosmological redshift. 

Our estimated SNRs suggest that coalescing black holes might well be the first sources detected 
by the LIGO/VIRGO network of ground-based interferometers. Because of their larger masses, they 
can be seen out to much greater distances (up to ~ 250 Mpc for M < 50Mq for the initial LIGO 
interferometers) than coalescing neutron star binaries (heretofore regarded as the “bread and butter” 
workhorse source for LIGO/VIRGO, visible to ~ 30Mpc by the initial LIGO interferometers). 

Low-mass BBHs (M < 30Mq for the first LIGO interferometers; M < 80Mq for the advanced; 
( 1 -|- 2 )M ^3x 10® M 0 for LISA) are best searched for via their well-understood inspiral waves; 
more massive BBHs must be searched for via their more poorly understood merger waves and/or 
their well-understood ringdown waves. A search for low-mass BBHs based on the inspiral waves, 
at a sensitivity level roughly half way between the hrst LIGO interferometers and the advanced 
LIGO interferometers, should be capable of finding BBHs out to ~ 1 Gpc. A search for massive 
BBHs based on the ringdown waves can be performed using the method of matched filters. If 
one wants the reduction in the event rate due to the discreteness of the template family to be no 
more than 10 %, then the number of independent templates needed in the search will only be about 
6000 or less. Such a search with the first LIGO interferometers should be capable of hnding BBHs 
in the mass range from about IOOM 0 to about 7 OOM 0 out to ~ 200 Mpc; with advanced LIGO 
interferometers, from about 2 OOM 0 to about 3 OOOM 0 out to 2 ~ 1; and with LISA, BBHs with 
10®M 0 (1 -|- z)M 3 X 10®M 0 should be visible out to 2 > 100. The effectiveness of a search 
based on the merger waves will depend on how much one has learned about the merger waveforms 
from numerical relativity simulations. With only a knowledge of the merger waves’ range of frequency 
bands and range of temporal durations, a search based on merger waves can be performed using a 
nonlinear hltering search algorithm. Such a search should increase the number of discovered BBHs 
by a factor of roughly 10 over those found from the inspiral and ringdown waves. On the other hand, 
a full set of merger templates based on numerical relativity simulations could further increase the 
number of discovered BBHs by a additional factor of up to about 4. 

RAGS numbers: 04.80.Nn, 04.25.Dm, 04.30.Db, 95.55.Ym 
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I. INTRODUCTION AND SUMMARY 

A. Coalescences of black hole binaries 

It has long been recognized that coalescences of bi¬ 
nary systems of two black holes could be an important 
source of gravitational waves both for the ground 

based interferometric detectors LIGO ^ and VIRGO Q] 
currently under construction, and also for the possible 
future space-based interferometer LISA The or¬ 

bits of binary black holes (BBHs) gradually decay from 
energy and angular momentum loss to gravitational radi¬ 
ation. Eventually, the holes coalesce to form a final black 
hole. For gravitational radiation reaction to successfully 
drive the binary to merge in less than a Hubble time, 
the initial orbital period must be < 0.3 days (M/M©)®/®, 
where M is the total mass of the binary; thus the critical 
orbital period is of order days for solar mass black holes, 
and of order years to hundreds of years for supermassive 
black holes (10®M© < M < lO^M©). 

The process of coalescence can be divided up into three 
more or less distinct phases: 

• An adiabatic inspiral, during which the gravita¬ 
tional radiation reaction timescale is much longer 
than the orbital period. The inspiral ends when 
the binary orbit becomes relativistically dynami¬ 
cally unstable at an orbital sraaration of r ^ 6M 
(in units where G = c = 1) |g,^. The gravitational 
waves from the inspiral carry encoded within them 
the masses and spins of the two black holes, some of 
the orbital elements of the binary and the distance 
to the binary 10 - 

• Towards the end of the inspiral, the black holes 
encounter the dynamical instability and make a 
gradual transition from a radiation-reaction driven 
inspiral to a freely-falling plunge |^,|ll|,^, after 
which, even if the radiation reaction could be 
turned off, the black holes would still merge. We 
will call the subsequent plunge and violent collision 
the merger phase. Gravitational waves from the 
merger could be rich with information about the 
dynamics of relativistic gravity in a highly nonlin¬ 
ear, highly dynamic regime about which we have a 
poor theoretical understanding today. 

• As the final black hole settles down to a stationary 
Kerr state, the nonlinear dynamics of the merger 
gradually become more and more describable as 
oscillations of this final black hole’s quasinormal 
modes ||^,Q . The corresponding emitted gravita¬ 
tional waves consist of a superposition of exponen¬ 
tially damped sinusoids. We will call the phase of 
the coalescence for which the emitted gravitational- 
wave signal is dominated by the strongest I = m = 
2 quasinormal mode signal the ringdown phase. 
The waves from the ringdown carry information 


about the mass and spin of the final black hole 
[pA 1T|. (Note that, for want of a better termi¬ 
nology, throughout this paper we consistently use 
coalescence to refer to the entire process of inspi¬ 
ral, merger and ringdown, and reserve the word 
merger for the phase intermediate between inspiral 
and ringdown.) 


In this paper we focus primarily on BBHs in which 
the masses of the two black holes are approximately the 
same, although we do also consider sources where one 
BH is much smaller than the other. We consider three 
different classes of BBHs: 


(i) Solar mass black hole binaries, i.e., binaries that 
are formed from massive main-sequence progenitor bi¬ 
nary stellar systems. Such BBHs are expected to have 
total masses in the range lOM© < M < 50M©, but not 
much larger than this. The rate of coalescence of solar- 
mass BBHs in the Universe is not very well known. How¬ 
ever, theory suggests that most of the BBH progenitor 
systems may not disrupt during the stellar collapses that 
produce the black holes, so that their coalescence rate 
could be about the same as the birth rate for their pro¬ 
genitors, about 1/100, 000 years in our Galaxy, or several 
per year within a distance of 200 Mpc 0 , ^8100]^ . 
Note that this coalescence rate is roughly the same as the 
expected event rate for what has traditionally been re¬ 
garded as the most promising source for ground based in¬ 
terferometers, coalescences of neutron star—neutron star 
(NS-NS) binaries ||01. The expected rate of NS-NS coa¬ 
lescences is more firmly known, however, since it is based 
on ex trapol ations from detected progenitor NS-NS sys¬ 
tems 01810. 

(ii) Intermediate mass black hole binaries, with total 
masses in the range 50Mq < M < (a few) x 10®M©. In 
contrast to the cases of solar mass black holes and su¬ 
permassive black holes (discussed below), there is little 
direct observational evidence for the existence of black 
holes in this mass range. Although there have been sug¬ 
gestions that the globular cluster M15 harbors a black 
hole of mass ~ 10®M© |0|, theoretical modeling com¬ 
bined with the most recent HST observations neither 
confirm nor rule out this possibility | 0 |. Despite the 
lack of observational evidence, it is plausible that black 
holes in this mass range are formed in the cores of glob¬ 
ular clusters, or in galactic nuclei in the process of for¬ 
mation of a supermassive black hole |0| . Simulations 
by Quinlan and Shapiro suggest that black holes with 
M ~ lOOM© — lOOOM© could be formed in the evolution 
of dense stellar clusters of main sequence stars in galac¬ 
tic nuclei 0], and that coalescences of binaries of such 
black holes could be possible en route to the formation 
of a supermassive black hole. 

We are prompted to consider these intermediate mass 
BBHs by the following consideration: even if the coa¬ 
lescence rate of intermediate mass BBHs is ^ 10“^ that 
of NS-NS binaries (which is thought to be ~ 10“®yr“^ 
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in our Galaxy as discussed above), these sources would 
still be seen more often than NS-NS binaries by LIGO’s 
initial and advanced interferometers, and thu s could be 
the first detected type of source. (See Sec. IE for further 
details.) 

(iii) Supermassive black holes binaries: There is a vari¬ 
ety of strong circumstantial evidence that supermassive 
black holes (SMBHs) in the mass range 10® — IO^Mq are 
present in quasars and active galactic nuclei, and also 
that ^ 25% - 50% of nearby massive spiral and elliptical 
galaxies harbor quiescent SMBHs See Ref. |^] 

for a review of this evidence. One of the main scien¬ 
tific goals of the LISA project is to detect and monitor 
various processes involving SMBHs, such as the capture 
of compact stars [pi[7| jT^J?7[ , ^ , and their formation . 
In particular, the coalescences of SMBH binaries that 
are formed in Galaxy mergers, in which the individual 
SMBHs are driven together by dynamical friction and 
gas accretion until gravitational radiation reaction takes 
over 1^ , have often been suggested a s a promis ing source 
for space-based interferometers |||,r 
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^ . Such co¬ 
alescences would be detectable throughout the observable 
universe with large signal to noise ratios 00- There is 
also observational evidence for SMBH binaries: wiggles 
in the radio jet of QSO 1928-1-738 have been attributed to 
the orbital motion of a SMBH binary as have time 
variations in quasar luminosities and in emission line 
redshifts The overall event rate is uncertain, but 

could be large (> 1/yr), especially if the hierarchical sce¬ 
nario for structure formation is correct 13 111 . 


B. Status of theoretical calculations of the 
gravitational-wave signal 

Detailed theoretical understanding and predictions of 
the gravitational waveforms /i+(t) and hx{t) produced 
in BBH coalescences will facilitate both the detection of 
the gravitational-wave signal, and the extraction of its 
useful information. In situations where a complete fam¬ 
ily of theoretical template waveforms is available, it will 
be possible to use the procedure of Wiener optimal fil¬ 
tering to search the interferometer data streams and to 
detect the gravitational-wave signal [0,|5| . The resulting 
signal-to-noise ratios (SNRs) can be larger than those ob¬ 
tainable without theoretical templates by a substantial 
factor; see Sec. Thus, while it is possible to detect the 
various phases of BBH coalescences without theoretical 
templates, such templates can greatly increase the effec¬ 
tive range of the interferometers and the event detection 
rate. (Accurate theoretical templates are also essential 
for extracting the maximum amount of information from 
detected signals |Q.) 

Such theoretical template waveforms are available for 
the inspiral and ringdown phases of the coalescence, but 
not yet for the merger phase, as we now discuss. 


For the inspiral phase of the coalescence, the gravi¬ 
tational waves and orbital evolution can be described 
reasonably well using the post-Newtonian approximation 
to general relativity. To date, inspiral waveforms have 
been calculated up to post-2.5-Newtonian order [^, and 
the prospects look good for obtaining waveforms up to 


post-3.5-Newtonian order 


Post-Newtonian tem¬ 


plates will be fairly accurate over most of the inspiral, 
the most important error being a cumulative phase lag 
[^,Q. This cumulative phase lag will not be impor¬ 
tant for searches for inspiral waves; template phasing 
error will be largely compensated for by systematic er¬ 
rors in best-fit values of the binary’s parameters, and 
the signals will still be found 42 By con¬ 


trast, template inaccuracies will be significant when one 
attempts to extract from the data the binary’s param¬ 
eters. In particular, post-Newtonian templates’ errors 
start to become very simificant around an orbital sep¬ 
aration of r ^ 12M [Q, well before the end of the 
inspiral at the dynamical orbital instability (r ^ 6M). 
Templates for the phase of the inspiral between roughly 
12M and 6M will most likely have to be calculated using 
methods other than the post-Newtonian approximation. 
Moreover, the methods of full blown numerical relativ¬ 
ity cannot be applied to this “Intermediate Binary Black 
Hole” (IBBH) phase, since the total time taken to evolve 
from 12M to 6M is about 1500M, too long for super¬ 
computer simulations to even contemplate evolving. Al¬ 
ternative analytic and numerical methods for calculating 
gravitational waveforms from the IBBH portion of BBH 
inspirals, based on the adiabatic approximation, are un¬ 
der development [^; it is likely that such alternative 
methods will be successfully developed and implemented 
before the gravitational-wave detectors begin their mea¬ 
surements 1 ^ . 

Waveforms from the dynamic, complicated merger 
phase of the coalescence can only be obtained from nu¬ 
merical relativity. Unlike mergers of neutron star bina¬ 
ries, BBH mergers are particularly clean in the sense that 
there is no microphysics or hydrodynamics to complicate 
simulations of the evolution, and external perturbations 
are negligible: the entire merger can be described as a 
solution to the vacuum Einstein equation |^. Finding 
that solution is not a particularly easy task: a major com¬ 
putational effort to evolve the vacuum Einstein equation 
for BBH mergers using massive computational resources 
is currently underway, funded by the National Science 
Foundation’s Grand Challenge program |0 , ^ . 

The ringdown phase of the coalescence can be accu¬ 
rately described using perturbation theory on the Kerr 
spacetime background . The gravitational waveforms 
from this phase are well understood, being just exponen¬ 
tially damped sinusoids. Thus, Wiener optimal filtering 
is feasible for searches for ringdown waves. 
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C. Purpose of this paper 


The principal purpose of this paper is to estimate, in 
more detail than has been done previously, the prospects 
for measuring gravitational waves from the three differ¬ 
ent phases of coalescence events (inspiral, merger and 
ringdown), for various different detectors, and for a wide 
range of BBH masses. We estimate in each case the dis¬ 
tances to which the different types of source can be seen 
by calculating expected SNRs. In particular, we deter¬ 
mine for each BBH mass and each detector whether a 
coalescence event is most effectively detected by search¬ 
ing for the inspiral portion of the signal, or the merger 
portion, or the ringdown portion. We also determine how 
much the availability of theoretical template waveforms 
for the merger phase could increase the event detection 
rate. 

Previous estimates of SNRs for ground-based interfer¬ 
ometers have focused on the inspiral |l|j4^ and ringdown 
1^,0 phases, and also focused on solar-mass BBHs. For 
space-based interferometers, previous estimates of SNRs 
from the merger phase were restricted to specific 

mass values and did not consider the ringdown portion 
of the signal. 

In a companion paper, we discuss in detail the use¬ 
ful information carried by the three phases of the 
gravitational-wave signal, and methods and prospects for 
extracting this information both with and without tem¬ 
plates for the merger phase |Q. 

In the following subsection we describe our calcula¬ 
tions and summarize the assumptions und erlying our es¬ 
timated signal-to-noise ratios. In Sec. IE we summarize 
our main results and con clusions, and an outline of the 
paper is given in Sec. IF. 


D. Estimating the signal-to-noise ratios: method 
and assumptions 


We calculate SNRs for three different types of inter¬ 
ferometer: initial and advanced ground-based interfer¬ 
ometers (LIGO/VIRGO), and the proposed space-based 
interferometer LISA. The noise spectra of the initial 
and advanced ground-based interferometers we took from 
Ref. dl , and that for LISA we obtained from Ref. @] . Our 
appro xim ate vers ions of these noise spectra are given in 
Eqs. (4.1) - (^), and are illustrated in Figs. |^ - || in 
Sec. 


We consider the following three different signal- 
detection methods: 


p, obtained from matched filtering is related to the grav¬ 
itational waveform h{t) measured by the interferometer 
and to the spectral density Sh{f) of the strain noise in 
the interferometer via |]^ 


p2=4 


\kf)? 

Shif) 


■df. 


(I.I) 


whe re h (f) is the F ouri er transform of h{t) defined by 
Eq. ( [2.31) . The SNR ( |l.lD depends, through the waveform 
h{t), on the orientation and position of the source rela¬ 
tive to the interferometer. In Sec. IIC we show that if we 


perform an rms average over source orientations and po¬ 
sitions (at a fixed distance), the rms SNR thus obtained 
depends only on the energy spectrum dE/df carried off 
from the source by the gravitational waves. The result¬ 
ing relationship between the waves’ energy spectrum and 
the rms angle-averaged SNR forms the basis for most of 
our calculations. It is given by [cf. Eq. (2.34)] 


{P^) = 


2(I-f z)2 


df- 


dE 


Jo ^ pSh{f) df 


[(1 + ^)/], ( 1 - 2 ) 


where z is the source’s cosmological redshift and D{z) is 
the luminosity distance from the source. 

In order for a signal to be detected, the waves’ mea¬ 
sured SNR must be larger than the detection threshold 


Pthreshold ~ y 2 ln(A^tart —times ) -|- 2 ln(A4hapes); (1-3) 

see, for example. Ref. [Q and also Sec. |IIB| . Here 
A4tart-times IS the number of independent starting times 
of the gravitational wave signal that are searched for in 
the data set, determined by the total duration of the 
data set (of order one year) and the sampling time. The 
quantity A4hapes is the total number of statistically inde¬ 
pendent waveform shapes in the set of signals that one is 
searching for . 

(ii) Band-pass filtering searches: For the merger phase, 
a complete set of theoretical templates may not be avail¬ 
able. It therefore may not be possible to use matched 
filtering, and other search methods will need to be em¬ 
ployed. Band-pass filtering, followed by setting a detec¬ 
tion threshold in the time domain, is a simple method of 
searching an interferomet er da ta stream for bursts of un¬ 
known form 1^ . In Sec. |H A| we derive an approximate 
relation between the S NR obtainable from band-pass fil¬ 
tering, and the SNR (O) obtainable from matched fil¬ 
tering, for any burst of waves, namely 


(i) Matched filtering searches: For those phases of the 
coalescence for which a complete set of theoretical tem¬ 
plates will be available (the inspiral, the ringdown, and 
possibly the merger), the method of matched filtering 
or optimal filtering can be used to search for the waves 
|p], ^5|j5^j5^j55|| . The theory of matched filtering is briefly 
sketched in Sec. HA, For any source of waves, the SNR, 


^ J band—pass V2TA7 \ 

/ optimal 


s 




(1,4) 


Here T is the duration of the burst an d A f is the band¬ 
width of the band-pass filter [cf. Eq. ( 2.14 )]. The quan¬ 
tity 2TAf is the dimension of the linear space of sig¬ 
nals being searched for, which is roughly the same as 
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the “number of cycles” of the g ravitational waveform. 
We use this formula in Sec. VIB to estimate band-pass- 
filter search SNRs for the merger waves from BBH co¬ 
alescences, by inserting on the right hand side the rms 
angle-averaged matched-filter search SNR (|l.2| ), and by 
making estimates of T and A/. 

(iii) Noise-monitoring, nonlinear filteri ng searehes: 
The traditional view has been that the SNR (IT) is about 
the best that can be achieved in the absence of theoret¬ 
ical templates; that is, that the gain in SNR obtainable 
from matched filtering is approximately the square root 
of the number of cycles in the gravitational wave signal. 
This view is based on the the assumption that the search 
method used in the absence of templates is band-pass 
filtering or som ething very similar. However, we sug¬ 
gest in Sec. 11B an alternative search method, motivated 
by Bayesian analyses and incorporating nonlinear filter¬ 
ing, which performs much better than band-pass filtering 
and in some cases almost as well as matched filtering. In 
essence, one monitors the noise level in the interferom¬ 
eter in a certain frequency band, over short timescales, 
and looks for statistically significant changes. The noise 
level is estimated by calculating the quantity 


the merger phase we make an educated guess of dE/df. 
Sec. Ill gives more details. 


1 

- / drs(^-hr)^ 
1 J-T/2 


(1.5) 


where T is the maximum expected duration of the signal, 
and s{t) is a suitably pre-filtered version of the interfer¬ 
ometer data stream. 

The efficiency of this noise-monitoring search method 
cannot usefully be described in terms of a signal-to-noise 
ratio, since the detection statistic is very non-Gaussian. 
Instead, its efficiency can be described in the following 
way. Let p denote the SNR that w ould be obtained if 
matched filtering were possible [Eq. 0)]- We use p as 
a convenient parameterization of the signal strength; as 
such, it is meaningful even in situations where matched 
filtering cannot be carried out. Then, a signal will be 
detected with high confidence using the noise-monitoring 
technique whenever p is larger than the threshold p*, 
where to a good approximation p* satisfies 

Ph< ~ 2 ln(A4tart —times ) -f Affainsln (l -f P*/A/'bins) ■ (1-6) 


Here A/bins = 2TA/ is the dimension of the signal space 
discussed above (or the number of independent frequ ency 
bins in the Fourier domain). The derivation of Eq. (l.f) 
is given in Sec. [11 B|. 

The relation ( |1.2| ) forms the basis of our SNR calcu¬ 
lations. We use the SNR thresholds 0 and 0) to 
deduce from the SNR values the detectability of the var¬ 
ious parts of the gravitational-wave signal. To calculate 
the SNRs, we also need to specify the waves’ energy spec¬ 
tra for the three different phases of the coalescence. As 
we now outline, for the inspiral and ringdown phases the 
waves’ energy spectrum is essentially known, while for 


Inspiral energy spectrum: We use the leading or¬ 
der expression for dE/df obtained using Newtonian 
gravit y su pplemented with the quadrupole formula |^] 
[Eq. ( 3.2C| )]. Strictly speaking, this spectrum describes 
the SNR that would be achieved by searching for New¬ 
tonian, quadrupole waves using Newtonian, quadrupole 
templates. The actual SNR obtained when searching for 
a real, general-relativistic inspiral waveform using post- 
Newtonian templates should deviate from this by only a 
few tens of percent |^. We terminate the spectrum at 
the frequency /merge = 0.02/M which is (roughly) the 
frequency of quadrupole waves emitted at the orbital dy¬ 
namical instability at r ~ 6M Q. For LISA, we assume 
that the measurement process lasts at most one year, 
and choose the frequency at which the inspiral spectrum 
starts accordingly. 

Ri ngdow n energy spectrum: The spectrum that we use 
[Eq. (3.24)] is determined (up to its overall amplitude) by 
the characteristics of the I = m = 2 quasi-normal ringing 
(QNR) mode of the final Kerr black hole. This mode is 
the most slowly damped of all QNR modes, so we expect 
it to dominate the last stages of the gravitational-wave 
emission. The QNR spectrum depends on three parame¬ 
ters: the quasi-normal modes’ frequency /qnr and damp¬ 
ing time T, and the overall amplitude of the quasinormal 
mode signal. Equivalently, the three parameters can be 
taken to be the mass M and dimensionless spin parame¬ 
ter a of the final black hole (which determine /qnr and r) 
and the total amount of energy radiated in the ringdown 
(which determines the overall amplitude). The spectrum 
is peaked at / = /qnr with width A/ ~ l/r. 

In our analyses, we (somewhat arbitrarily) assume that 
a = 0.98. It seems likely that in many coalescences the 
spin of the final black hole will be close to maximal, since 
the total angular momentum of the binary at the end of 
the inspiral is ~ 0.9M^ when the individual black holes 
are non-spinning [^ , and the individual black hole spins 
can add to this. Exactly how close to maximally spinning 
the final black hole will be is a matter that probably 
will not be decided until supercomputer simulations—or 
observations—settle the issue. In any case, the ringdown 
SNR values that we obtain depend only weakly on our 
assumed value of a [cf. Eq. ( |A14 )], for fixed total energy 
radiated in the ringdown. 

The overall amplitude of the ringdown signal depends 
upon one’s delineation of where “merger” ends and “ring- 
down” begins, which is somewhat arbitrary. For equal- 
mass BBHs, we assume a value of the overall ampli¬ 
tude that corresponds to a total radiated energy in the 
ringdown of 0.03M—a radiation efficiency of 3%. This 
number is based on a back-of-the-envelope, quadrupole- 
formula-based estimate of the QNR mode’s amplitude 
when the distortion of the h orizon of the black hole is 
of order unity {cf. Sec. HID). Although this radiation 


efficiency seems rather high, there have been numerical 
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simulations of the evolution of distorted, spinning black 
holes in which the ringdown waves carry away ^ 3% of 
the black hole’s total mass energy [Q. 

For non equal-mass black holes, we assume that the to¬ 
tal energy radiated is F(/r/M)0.03M, where F{fj,/M) = 
(4/r/M)^ and /r is the reduced mass of the binary. This 
function gives the correct results for the equal-mass case 
and also gives the correct scaling law in the regime 
/i M; for general mass ratios the scaling law is proba¬ 
bly a good approximation. 


Merger energy speetrum: Realistic merger energy spec¬ 
tra will vary substantially from event to event (depend¬ 
ing on the initial spins of the inspiraling black holes). 
Currently, we have very little concrete information about 
such spectra, pen ding su percomputer simulations of BBH 
mergers. In Sec. niB| we describe various circumstan¬ 
tial pieces of evidence, culled from the literature, rele¬ 
vant to merger spectrum. Based on that evidence, we 
adopt the following crude model for equal-mass BBHs: 
we assume a flat spectrum dE/df = const extending from 
the frequency /merge = 0.02/M of quadrupole waves at 
the end of inspiral to the quasinormal ringing frequency 
/qnr = 0.13/M, with amplitude such that the total radi¬ 
ated energy in the me rger i s 10 % of the total ma ss ene rgy 
of the spacetime [Eq. ( 3.14 )]. We outline in Sec. IIIB two 
different “handwaving” arguments which suggest that in 
favorable cases the merger radiation efficiency may be as 
high as our assumed value of ~ 10%. One of these ar¬ 
guments, due originally to Smarr and explored by 
Detweiler , is based on extrapolation of perturbation 
theory results; the other argument is based on angular 
momentum conservation. 

Our assumed radiation efficiencies of 3% and 10% for 
the ringdown and merger phases should perhaps be inter¬ 
preted as reasonable upper bounds that could be achieved 
in favorable cases, rather than as best-guess estimates. 
We note that numerical simulations that have been per¬ 
formed to date (which are restricted to axisymmetric sit¬ 
uations) generally yield lower radiation efficiencies than 
we have assumed p™; moreover, these axisymmetric sim¬ 
ulations generally find that ringdown waves carry most 
of the radiated energy. In Sec. IIIB we argue that the 
radiated energy in the merger phase could be boosted by 
the lack of symmetry in generic black hole mergers, and 
especially by the individual black holes’ spins (if these 
spins are large). 

For non-equal mass BBHs, we again simply reduce 
the energy spectrum by the factor F{fj,/M) = (An/M)'^, 


while the upper and lower frequencies /merge 
taken to be independent of /r. 


and /q 


are 


E. Signal-to-noise ratios: results and implications 


By inser ting our assumed energy spectra ( 3.20| ), (3.14) 
and (3.24) into Eq. (0), we obtain optimal-filtering 
SNRs for the three different phases of BBH coalescences 


as a function of the redshifted total mass (1 -|- 2 )M of the 
binary. The results are shown in Figs. || and ^ and 
formulae summa rizing the results are given in Appendix 
A. Also in Sec. HIE we estimate that, for the merger 
waves, the number of independent frequency bins A/bms 
which characterize the signal falls roughly in the range 
10 < A/bins ^ 30, and that a conservativ e upp er-bound is 
~ 60. We use this upper bound in Sec. VIB to estimate 
the SNR threshold (0) for merger waves using noise¬ 
monitoring searches when templates are unavailable. We 
discuss the implications of these SNRs and SNR thresh¬ 
olds in Sec. VI; here we summarize our main conclusions: 


• Ground-based interferometers can study black- 
hole mergers in the mass range (a few)M 0 to ~ 
2 OOOM 0 ; LISA, by contrast can study mergers in 
the mass range 10^Mq ^ (1 -I- z)M < 10®M 0 . 


• Ground-based interferometers can do moderate 
SNR (a few tens), moderate accuracy studies of 
the dynamics of merging black holes. LISA, by con¬ 
trast, can do high SNR (a few x 10^), high-accuracy 
studies. 


• Goalescing black holes are fairly likely to be the 
first sources detected by the LIGO/VIRGO net¬ 
work of ground-based interferometers: because of 
their larger masses, they can be seen out to much 
greater distances than coalescing neutron star bi¬ 
naries. (With the initial LIGO interferometers, bi¬ 
nary black holes with M ^ 50Mq can be seen out 
to ~ 250 Mpc, whereas binary neutron stars can be 
seen out to ~ 30 Mpc p^). Estimates of the num¬ 
ber of BBH coalescences per NS-NS coalescence by 
experts in binary evolution range from 1/10 to 
1/300 1^1; the distance gain factor for BBHs could 
easily compensate for the reduced birth rate. 


• Low-mass BBHs [M // 30Mq for the first LIGO 
interferometers; M < 80Mq for the advanced; 
(1 -I- z)M < 3 X 1 O®M 0 for LISA] are best searched 
for via their well-understood inspiral waves; while 
more massive BBHs must be searched for via their 
more poorly understood merger waves and/or their 
well-understood ringdown waves. 


• When the interferometers’ sensitivities has reached 
a level roughly half way between the first LIGO 
interferometers and the advanced LIGO interfer¬ 
ometers 1 ^, a search based on the inspiral waves 
should be capable of finding low-mass BBHs out to 
about 1 Gpc. 


• A search for massive BBHs based on the ring- 
down waves can be perf ormed using matched fil¬ 
tering. We show in Sec. VIA that the number of 
independent templates needed for such a search is 
about 6000 or less, assuming that one wants the 
event rate reduction due to discreteness of the tem¬ 
plate family to be no more than 10%. With the 
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first LIGO interferometers, such a search should be 
capable of finding equal-mass BBHs in the mass 
range IOOMq to 700Mq out to about 200 Mpc; 
with advanced LIGO interferometers, BBHs with 
200Mq < M(1 + z) < 3OOOM0 should be de¬ 
tectable out to z ^ 1; and with LISA, BBHs with 
IO^Mq <{1+ z)M < 3 X 10®Mo should be visible 
out to z > 100. These distances are reduced by a 
factor ~ (4/i/M) for non-equal-mass BBHs. 

• The effectiveness of a search based on the merger 
waves will depend on how much one has learned 
about the merger waveforms from numerical rela¬ 
tivity simulations. With only a knowledge of the 
merger waves’ range of frequency bands and range 
of temporal durations, a search based on merger 
waves can be performed using the noise-monitoring 
search algorithm discussed above. Such a search 
should increase the number of discovered BBHs by 
a factor ~ 10 over those found from the inspiral 
and ringdown waves. On the other hand, a full set 
of merger templates based on numerical relativity 
simulations could further increase the number of 
discovered BBHs by an additional factor of up to 
- 4. 


F. Organization of this paper 


The organization of this paper is as follows. In Sec. || 
we discuss in detail the three methods of searching for 
gravitational wave signals referred to ab ove. In Sec. p C] 
we derive the general formula [Eq. ( 2.34 )] discussed above 
for the optimally-filtered SNR of a gravitational-wave 
source at a cosmological distance, in terms of the source’s 
emitted spectrum, redshift and luminosity distance, and 
the detector’s noise spectrum. This formula is one un¬ 
derpinning for the BBH SNR computations of Appendix 
A and Sec. In S ec. HA we derive the approximate 
formula [Eq. ( |2.I4 )] discussed above for the SNR ob¬ 
tainable by band-pass fil terin g in the absence of theoret¬ 
ical templates. In Sec. 1^ we describe our suggested 
noise-monitoring sea rch, and d erive the relevant detec¬ 
tion thres hold [Eqs. ( 2.30D and ( 2.31 )]. 

In Sec. Ill A we split the waveform from a coalescing 


equal-mass black hole binary into three successive epochs: 
inspiral, merger and coalescence. Then we estimate the 
frequency bands and spe ctra for eac h o f t hes e e poch s: 
the merger phase in Sec. tllB [Eqs. (|3.^ , (3.5), ( 3.14 )[; 
the inspiral phase in Sec. IHC [Eq. (T^)[; and the ring- 


down phase in Sec. IH 1 [Eqs.(|3.24| ) - (|3.26|) [. We have 


already summarized, above, the assumptions that under¬ 
lie our assumed merger and ringdown spectra, and the 
uncertainties in those ass umpt io ns; w e give a more de¬ 


tailed description in Secs. HIB - HID, A subtlety in the 


derivation of the ringdown energy spectrum is explained 


in Appendix B. In Sec. HIE we make rough estimates of 


the range of the bandwidths and temporal durations of 
merger waveforms. 

In Sec. 0 we devise a simple piece-wise power-law ana¬ 
lytic fit [Eq. Q] to the noise spectra of an initial LIGO 
interferometer, an advanced LIGO interferometer, and a 
space-based LISA interferometer. This single formula, 
by adjustment of its parameters, can describe all three 
interferometer types. 

In Appendix A, we insert the interferometer noise 
spectra from Sec. IV, and the estimated BBH inspiral, 


merg er, and ringdown spectra from Secs. 1^ , pT^ and 


HID, into the general SNR formula (2.34). That for¬ 


mula’s frequency integral is then performed to compute 
the angle-averaged, optimal-filtering SNR for each type 
of interferometer and each phase of BBH coalescence. 
The resulting SNRs for an initial LIGO interferometer, 
an advanced LIGO interferometer, and a LISA interfer¬ 
ometer, are expressed in terms of the binary’s total mass 
M, cosm ologi cal redshift z, and luminosity distance D. 

In Sec. VB| , we graph the SNR formulae derived in Ap¬ 
pendix A (Fig. Hlfor an initial LIGO interferometer; Fig. || 
for an advanced LIGO interferometer; Fig. for LIS A). 
We give intuitive insight into these SNRs in Sec. V_A by 
re-expressing the power SNR for a source as 


P^= d(ln/)[/lchar(/)//ln(/)]^ 


(1.7) 


where /ichar(/) is the source’s “characteristic amplitude” 
as a function of frequency, and hn{f) is the detector’s 
rms noise in a bandwidth equal to frequency for sources 
with random orientations. We give plots of /ichar(/) and 
hn{f) for five specific examples of binaries — a IOMq — 
IOMq binary at distance D = 200 Mpc and a 15Mq — 
15Mq binary at cosmological redshift z = 1 (Fig. |^); a 
50Mq — 50Mq binary at cosmological redshift z = 0.5 
(Fig. H); a 5 X 10®Mq — 5 x 10®M 0 binary at z = 5 and a 
2.5 X 1O^M0 — 2.5 X 1O^M0 binary at z = 1 observed by 
LISA (Fig. H). For each example, we explain, in terms of 
the plots, why the inspiral, mercer, and ringdown SNRs 


i I OT |. 

and H are the foundation 


have the values shown in Figs. 

The SNR graphs in Figs. | 
for our conclusions, summarized above, about what fea¬ 
tures of which binaries should be observable with which 
interferometers. A detailed discussion of these conclu¬ 


sions is given in Sec. VI. In Sec. VIA, we estimate the 
number of templates required for a search for ringdown 
waves based on matched filtering, estimate the SNR de¬ 
tection thresholds, and hence the range of the various in¬ 
terferometers for ringdown waves. In Sec. VIB we exam¬ 


ine the prospects for searches for BBHs via their merger 
waves, both with and without templates. 
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II. DERIVATION OF GENERAL FORMULAE 
FOR SIGNAL-TO-NOISE RATIOS AND 
DETEGTION THRESHOLDS 

In this section we discuss in detail the various signal- 
search methods which were brie fly described in the In¬ 
troduction. We star t in Sec. IIA by deriving the approx¬ 
imate relation ( |1.4D between the SNR achievable using 
matched filtering searches for signals and the SNR ob¬ 
tainable via band-pass filtering searches. This approxi¬ 
mate relation is closely related to the standard lore that 
the gain factor due to matched filtering is the square 
root of the number of cycles in the waveform 10,11, which 
strictly speaking is applicable only to waveforms that are 
alm ost m onochromatic. For completeness, we disc uss in 
Sec. [I A both the standard lore relation [Eq. jil] and 
the approximate more general relation [Eq. CT)T- 
In Sec. IIB we describe our proposed noise-monitoring 
search method, and derive the detection threshold ( |1.6| ) 
discussed in Sec. Fi nally , in Sec. I^ we derive 

the general formula [Eq. ( 2.34 )[ for the angle-averaged, 
optimal-filtering SNR for a gravitational-wave source, 
which was discussed in the Introduction. 


A. Searches for gravitational-wave bursts: band-pass 
filtering and matched filtering 

Consider the situation where 

some arbitrary gravitational-wave burst h{t) is present 
in the data stream s{t), so that 


s{t) = h{t) + n{t), 


( 2 . 1 ) 


where n(t) is the noise. If one integrates a filter K{t) 
against the data stream s(t) to produce a number, Y = 
f K(t)s(t) dt, then the standard definition of the signal 
to noise ratio is Q 

S expected value of Y when signal present 
JV rms value of Y when no signal present 

(y) 




4/“d/5i 

'kf)*Kif) 

\I^Io^df 

Kif) 

^ Shif) 


( 2 . 2 ) 


see, e.g., Refs. [ p3y54| ]. Here tildes denote Fourier trans¬ 
forms according to the convention 

/ OO 

(2.3) 

-OO 

and Shif) is the power spectral density of strain noise in 
the detector [Q. 

Now, consider searching for a signal h{t) when the 
only information about it that one has is its approxi¬ 
mate bandwidth in the frequency domain. Perhaps the 


simplest search algorithm one could use to search for the 
signal is to choose for K{t) the following band-pass filter: 


K(f) = e2"^‘--0(A//2 - 1/ - /eharl). 


(2.4) 


Here 0 is the step function and fstart is the starting time 
of the filter. This filter chops out all the data in the 
frequency domain except that in a bandwidth A/ about 
a characteristic central frequency /char @] • Suppose that 
the frequency interval has been chosen wisely, so that 
the signal h{t) has negligible power outside the interval. 
Then h{f) can be ta ken to va nish outside the chosen 
bandwidth, and Eqs. (2.2) and (|2.4D yield 


^(^start) 


band—pass filter 


/kfdfShif) 

/char h{t start) 


A/ h rms(/char) 


(2.5) 


where /irms(/) = yfShJj) is the rms fluctuation in the 
noise at frequency / in a bandwidth equal to /. The 
starting time of the filter fstart is then varied to give 
the maximum filter output Y, which is achieved at some 
value fbest of tstart- At t his maximum overlap time, the 
SNR is given by Eq. ( |2.5D with fstart replaced by fbest- In 
particular, for broadband signals for which A/ ~ /char 
(as opposed to Af «C /char), Eq. (^) simplifies to the 
standard result [M 




band—pass filter 


^rms (/char) 


( 2 . 6 ) 


By contrast, if the shape of the signal is known one can 
use the well known optimal filter K{f) = h(f)/ Shif), for 
which the SNR squared is, from Eq. (^) IsF 


'■■(I) 


= 4 


optimal filter 


f 


IM/)I^ 

Shif) 


■df. 


(2.7) 


A crucial element of both optimal filtering searches and 
most especially band-pass filtering searches with ground- 
based interferometers is the use of coincidencing between 
different interferometers to circumvent the effects of non- 
Gaussian noise bursts |35|. Coincidencing between the 4 
interferometers in the LIGO/VIRGO network (the Han¬ 
ford 2 km, Hanford 4 km, Livingston 4 km and Pisa 3 
km interferometers) should be sufficient to achieve this. 
To be conservative, our assumed detection thresholds for 


the SNR values (which we discuss in Sec. VI) are based 
on combining just the two LIGO 4 km interferometers, 
albeit with assumed Gaussian statistics. 

The standard lore appr oxim ate relationship between 
the optim al-fi ltering SNR (^) and the band-pass filter¬ 
ing SNR (2.6) can be obtained as follows. Consider the 
special case of a waveform that is quasi-monochromatic, 
i.e., of the form 































hit) = 

^ainp (t) cos [$(t)] , 
amplitude /lamp (/) 


( 2 . 8 ) 


where the amplitude /iamp(/) and instantaneous fre¬ 
quency [given by 27r/(/) = d^/dt] are slowly evolving. 
Using Eq. (2.7) and the stationary ph ase a pproximation 
to the Fourier transform of the signal (|2.8|) yields 


= J ricycif) 


mr 


sify 


(2.9) 


where ricydf) — P/f is the number of cycles spent 
within a bandwidth A/ ~ / centered on /, and /(/) 
is the time at which th e gr avitat iona l-wave frequency is 
/. By comparing Eqs. (2.6) and (2.9) it can be seen that 


ncydf) is the gain factor in SNR squared for optimal fil¬ 
tering over band-pass filtering, per logarithmic interval 
in frequency [||. 

For general sign als /i(/) which are not quasi- 
monochromatic, Eq. (|2.8| ) does not apply. However, we 
can derive an approximate formula for the SNR ( |2.6| ) for 
general signals as follows. Approximating Shif) to be 
constant in Eq. (2.7) gives 


S' 


V _ 2 7°“ 

^ J optimal filter S/j(/char) J —c 
« S/charT 


dt [hit)Y 




/irms (y*char)^ 


( 2 . 10 ) 


where h is a rms average of h{t) and T is th e effe ctive 
duration of the signal. Comparing Eqs. (2.10) and (2.1;) 
we find that 


information or assumptions about the signal and not the 
signal itself. This is because for band-pass filtering we 
must choose at the start some bandwidth A/. Since the 
true bandwidth of the signal will vary from one signal 
to another and will not be known in advance, the true 
bandwidth of the signal will in general be somewhat less 
than A/ if we have chosen A/ wisely. The quantity A/bins 
depends not on the true bandwidth of the signal but on 
our assumed bandwidth A/. 


The first factor on the right hand side of Eq. ( 2.11 ) 
is the ratio between the peak strain amplitude /i(/best) 
in the time domain and an rms value h of this strain 
amplitude. By defining the effective duration T of the 
signal to be given by 


dt [hit)f = T hithesp 


(2.13) 


this facto r red uces to unit y. With this interpretation of 
T in Eq. ( 2.12 ), Eq. ( ^.11 ) reduces to the formula 


iS/N) 


band—pass filter 


1 


(2.14) 


(*^/-^)optimal filter V-MDins 

discussed in the Introduction. We use this formula in 


Sec. VI 


The SNR threshold ^threshold appropriate for matched 
filtering is approximately given by 


erfc(/?threshold/'N/2) — 


A4tart —times A/shapes 

which to a good approximation reduces to 


(2.15) 


(‘^/A^)band_pass filter /l(/best) 1 

(^/^)optimal filter ^ V"A/bins 

where the “number of a priori frequency bins” 

Afbins = 2TA/. 

The reason for this terminology is as follows, 
that we are trying to detect a signal whose total dura¬ 
tion we know to be less than or equal to T, and whose 
spectrum we know to lie in some frequency interval of 
bandwidth A/. When we discard the measured data out¬ 
side of these time and frequency intervals, the remaining 
relevant data is described by A/bins = 2TA/ real Fourier 
coefficients (or equivalently, A/bms frequency bins). The 
quantity A/bins is the number of independent real vari¬ 
ables that parameterize the space of signals for which we 
will search. 

This notion of number of a priori frequency bins is 
closely related to the notion of number of cycles in the 
waveform. The number of waveform cycles will be ap¬ 
proximately Afcyc ~ T/char, and so for a broadband burst 
for which A/ « /char we have N'cyc ~ A/bins- There is a 
distinction between the two concepts, however: the num¬ 
ber of cycles is intrinsic to the signal, whereas the num¬ 
ber of frequency bins in part characterizes our a priori 


( 2 . 11 ) 


IS 


( 2 . 12 ) 


Simnnsp 


^threshold ~ \J 2 lll(A4tart —times /c) 2 ln(A/shap es) 5 

(2.16) 

[c/. Eq. (^/^) where the factor of e, which does not 
have a large effect, was omitted]. Here the number 
A/shapes “ A^shapes(/^threshold) 1® the number of Statisti¬ 
cally independent waveforms with SNR < pthreshoid in 
the set of signals to be searched for [^; Eq. ( 2.16 ) must 
be solved self-consistently to determine pthreshoid- 


B. Searches for gravitational-wave bursts: noise 
monitoring 

In this section we describe our suggested noise¬ 
monitoring method for searching for gravitational wave 
bursts of unknown form; more details can be found in 
Ref. iQ. In essence, the method consists of monitoring 
over short timescales the total rms noise in the detector 
output in the frequency band in which the signal is ex¬ 
pected, and waiting for statistically significant changes 
in one’s estimate of the noise power. 

Suppose that the maximum expected signal duration 
is T, and that the interferometer output is s(/). We con¬ 
struct a quantity Qit) in the following way. First, focus 
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attention on the data stream s(t) in the time interval 
t — T/2<T<t + T/2. Because the data stream is in 
fact discrete and not continuous, this interval of data can 
be represented by the numbers 

Sj = s{t-T/2 + jAt) (2.17) 

for 0 < j < fVtotai = T/At, where At is the sampling 
time. From Eq. ( ^.l[ ) we have 

sj = hj+nj, (2.18) 


Thus, the quantity is the optimal filtering signal-to- 
noise squared per data bin. Note that throughout this 
subsection, we use p as a convenient parameterization of 
the signal strength, which is meaningful even in situa¬ 
tions where templates are not available and where op¬ 
timal filtering cannot be carried out. Band-pass filter¬ 
ing (of a pre-whitened data stream) approximately cor¬ 
responds in this language to calculating the statistic 

PBP = max—. (2.24) 

j aj 


where hj is the gravitational wave signal and rij is the 
noise. Now because the interferometer noise is colored, 
the noise matrix 




(2.19) 


will not be diagonal. Here, angle brackets denote en¬ 
semble averaging over realizations of the noise. If one 
performs an FFT just of this finite stretch of data, the 
noise matrix on the new basis will not be diagonal ei¬ 
ther because of aliasing effect s. However, it is possible to 
diagonalize the matrix ( 2.19| ) and change to a basis on 
which the noise is diagonal. We will denote this new ba¬ 
sis by capital Roman letters I, J, K. The data points s/ 
on this new basis can be chosen to correspond approxi¬ 
mately to frequ encies // = I/T, 1 = 1, —1, 2, —2,... [ ^ . 
Equation (2.19) can now be replaced by 


{mnj) = Sija], 


( 2 . 20 ) 


The data sj extend up to some high frequency (of or¬ 
der several kHz) determined by the sampling time. We 
next discard all data above some upper cutoff frequency. 
(Note that we have effectively performed band-pass fil¬ 
tering of the data, since the restriction to a segment of 
length T in the time domain removes frequency compo¬ 
nents at / < 1/T.) Let us denote by A/bins the total 
number of data points remaining, which will be approxi¬ 
mately given by 


Nuns = 2TAf 


( 2 . 21 ) 


where A/ is the bandwidth of our effective band-pass 
filter. 

In terms of this basis, matched filtering consists of cal¬ 
culating, for each trial waveform shape hj, the quantity 


Y.jSjhj/a'] 

vTT^- 


( 2 . 22 ) 


(We are assuming here that all the trial waveform shapes 
have duration less than T and in the frequency domain 
have most of their power within the bandwidth A/). We 
introduce the notation pi = hi/aj', then the optimal 
filtering SNR (H)is given by 


= 


Ep? = e3- 

I I ^ 


(2.23) 


This will hav e an expected value of roughly pj VA^bins 
[cf. Eq. (2.14)] if the signal is not sharply peaked in fre¬ 
quency but instead is spread out over the bandwidth A/. 

The statistic used to diagnose the presence of a signal 
is 


Q{t) — — Abins + 


J —■A^bins 


/2 ,2 

II. 
2 ■ 


J — Nbins / 2 


a 


(2.25) 


J 


This statistic is essentially an estimate of the apparent 
rms noise power in the given bandwidth over the given 
time interval, up to an additive constant. The additive 
constant, —A/bms, is chosen so that when no signal is 
present, {Q{t)) = 0 and thus Q{t) random walks over 
positive and negative values. When a signal is present, 
Q{t) will with high probability be large and positive. One 
monitors Q{t) as a function of time and sets a detection 
threshold such that when Q{t) exceeds the threshold, the 
probability that the data stream contains no signal is very 
low (see below). Note that this search method constitutes 
a type of nonlinear filtering. 


The technique of monitoring the statistic ( 2.25 ) is 
closely related to two commonly used techniques in radio 
astronomy. In the first such technique, observers sum the 
power from frequency bins which are expected to contain 
harmonics of the signal they are trying to detect. This 
procedure is not quite as good as coherently combining 
the signal from all the frequency bins but is much eas¬ 
ier, computationally. The second technique is use¬ 
ful when one is looking for periodic signals in a data 
train that is too long to Fourier transform. One splits 
the data train into some number N of shorter data seg¬ 
ments, takes the FFT of each shorter segment, and adds 
the FFTs incoherently (he.., adds the individual power 
spectra). Again, this is not the optimal search method, 
but it is often a useful thing to do given finite compu¬ 
tational resources. One difference between our suggested 
noise-monitoring search method and the radio astronomy 
techniques is the following: One a dds the frequency bins 
incoherently in the statistic ( ^.25 ) because the phase re¬ 
lationships are unknown; in the radio astronomy con¬ 
texts, one adds the frequency bins incoherently to save 
on computation time. But the techniques themselves are 
conceptually very similar. 

We now turn to a derivation of the efficiency and 
performance of our suggested noise-monitoring search 
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method. It is straightforward to show that when a signal 
is present 


(Qit)) = 


(2.26) 


and that 


([Q(i) - (Q(<))]") = V + 2Afbins. (2.27) 


When no signal is present, Eqs. ( ^.26 ) and ( ^.27 ) continue 
to hold with p = 0. These equations show that a signal 
should be detectable in the regime 




(2.28) 


as well as at l arger p, since in the regime ( 2.28] ) the ex¬ 
pected value ( 2.26 ) of Q is large compared to its rms 
value in the absence of a signal. By contrast, a signal 
is detectable using band-pass filtering only in the regime 
p ^ from Eq. ( ^.24 ) and associated discussion. 


The noise-monitoring method thus extends the domain 
of detectable signals to include the regime ( ^.28 ). 

T he a pproximate SNR threshold predicted by 
Eq. ( 2.27| ) is correct in order of magnitude, but to obtain 
an accurate SNR threshold it is necessary to calculate the 
full probability distribution for the statistic Q, since its 
statistical properties are very non-Gaussian. This prob- 
abilit y dis tribution is given by, from Eqs. (2.18), ( 2.2(\ ) 
and (2.25), 


P[Q{t) > Qo] = 


r[A/bins/2, (Qo + A/'bins)/2] 
r(Afbins/2) 


(2.29) 


where r(- • •, • • •) is the incomplete Gamma function and 
r(- • •) is the usual Gamma function. Suppose now that 
we examine a number A4tart-times of different starting 
times t, and that we w ish to hnd the number Qq such that 
the probability ( 2.29| ) of Q{t) exceeding Qq for any t, in 
the absence of a signal, is some small number e ^ 10“^. 
This threshold Qo is obtained by solving 


r[A/bins/2, (Qo + A/bins)/2] 
r(Afbins/2) 


A4tart — 


(2.30) 


start—times 


From Eq. o, this threshold wil l be exceeded by a 
signal whenever the signal strength (^.23) satisfies 


> P* = \/Qo- 


(2.31) 


Equations ( 2.30| ) and ( ^.3lD determine the threshold p* 
as a function of the parameter s e, A 4tart-times; and A^bins; 
we use these formulae in Sec. VIB. 


We now discuss an approximate formulae for the 
threshold p* which gives an insight into the nature of 
the nonlinear filtering algorithm. It can be shown |5|] 
that for A/bins 1, p* can be obtained to a good approx¬ 
imation by solving the equation 


— 2 ln(A/start—times/c) -f A^bins In (l 4- P^/J^hins) : 


(2.32) 


[cf. Eq. ( |l.6| ) discussed in the Intro duct ion]. It is possible 
to understand the origin of Eq. ( 2.32| ) in the following 
simple way. 

The total number A4hapes.max of waveform shapes of 
duration less than or equal to T, with bandwidth A/ 
and with SNR < p which are distinguishable in the in¬ 
terferometer noise is approximately given by ||36| 


ln(A4hap 


es.max 


) ~ In (l + P^/A/bins) ■ (2.33) 


Let us denote by Ad the manifold of all possible wave¬ 
forms satisfying these criteria — duration less than or 
equal to T, bandwidth A/ and SNR < p. Not all of 
the waveforms in this manifold correspond to physically 
interesting BBH merger waveforms; the set Admerge of 
BBH merger waveforms will form a sub-manifold of Ad. 
Consider now the hypothetical limit in which the sub¬ 
manifold Ad merge becomes so large (perhaps curving back 
and intersecting itself) that, when smeared out by the in¬ 
terferometer noise, it effectively fills the entire manifold 
Ad. In this limit, the number Afshapes of BBH merger 
waveform shapes for which we are searching approaches 
A4hapes,max- It is dear that in this limit, knowledge of 
the waveform shapes is not useful in searching for the 
waves, and that a search method which simply seeks the 
maximum overlap between the measured signal and all 
possible waveform shapes (thus not requiring templates) 
would perform just as well as matched filtering. This is 
essentially what the nonline ar se arch method does. Using 
this insight, the thr eshol d ( 2.32|) can be obtained simply 
by combining Eqs. ( ^T6l ) and (^.331). 

The above derivation is based on frequentist statistics. 
In Ref. a Bayesian analysis is outlined of the de¬ 
tection of gravitational wave signals of unknown form 
which automatically identifies the statistic Q(t) as opti¬ 
mal, and which also approximately reproduces the detec¬ 
tion threshold p*. 

In practice, this search method would be combined 
with coincidencing between interferometers to achieve 
high detection reliability and to reduce the effects of non- 
Gaussian noise, a s is t he case with band-pass filtering as 
discussed in Sec. IIA. Matched filtering could be more 


efficient than the noise-monitoring method at combat¬ 
ing non-Gaussian noise via coincidencing, for the follow¬ 
ing reason. When coincidencing with templates, one can 
demand that the SNR in each interferometer be above 
the appropriate threshold, and that the signal-parameter 
values deduced in each interferometer be consistent with 
each other. For the noise-monitoring searches, one can 
only demand that the SNR in each interferometer be 
above the appropriate threshold. Hence, matched fil¬ 
tering has more discriminating power against situations 
in which all the interferometers have moderately large 
non-Gaussian noise spikes somewhere in the relevant 
time window. Non-Gaussian noise may therefore make 
the less-discriminating noise-monitoring search perform 
somewhat worse in practice, relative to matc hed fi ltering 
searches, than is indicated by the threshold (2.32|). 
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C. Signal-to-noise ratio for matched filtering in 
terms of waves’ energy spectrum 


In this section we derive the relation [Eq. (2.34)] be¬ 
tween the expected value of the optimal-filtering SNR 
(2.7), and the energy spectrum of gravitation al w aves 
emitted by the source. In general, the SNR (H) for 
a burst of waves depends on the details of the gravita¬ 
tional waveform, on the orientation of the source with 
respect to the interferometer, and on the direction to 
the source. By contrast, the quantity (p^), the average 
of the squared SNR over all orientations of and direc¬ 
tions to the source, depends only on the total energy 
per unit frequency dE/df carried off from the source 
by the waves. Specifically, consider a gravitational-wave 
source located at a cosmological redshift 2 ; and corre¬ 
sponding luminosity distance D{z). Let the locally mea¬ 
sured frequency of the waves near the source be /e, re¬ 
lated to the frequency / measured at the interferometer 
by / = /e/(I + z). Let the locally measured energy spec¬ 
trum of the waves be dE^/dfe{fe)- Then the orientation- 
averaged SNR squared measured at the interferometer is 
given by 


{P^) = 


2{lEzf 


I 


df- 


1 


57T^D{zy Jo ■' pShif) df, 


^[(I 


-z)f]. (2.34) 


Note that the relation (2.34) refers to an angle- 
averaged SNR obtained from an rms average of sig¬ 
nal amplitudes over different possible orientations of the 
source and interferometer. This averaging convention dif¬ 
fers from that adopted in Refs. 00 - There, the angle- 
averaged SNR is defined to be a cube root of an average 
of cubed signal amplitudes, rather than an rms average. 
That “cube root of a mean cube” averaging method is the 
appropriate method for calculating the expected event 
detection rate Q. As a result, the SNR formulae used in 
Refs. are a factor of -^3/2 larger than the formu¬ 

lae used in this paper, the factor of -^3/2 being an ap¬ 
proximation to the effect of the different angle-averaging 
methods. 


We now turn to the derivation of Eq. ( 2.34 ). Consider 
first the case where the source of gravitational waves 
is sufficiently close that cosmological effects can be ne¬ 
glected. Let the source be at a distance r from the detec¬ 
tor and at a location (0, ip) on the sky. Let (i, (3) denote 
the direction towards the detector (spherical polar an¬ 
gles) with respect to a set of Cartesian axes centered at 
and determined by the source. Let the two independent 
polarizations of the strain amplitude at the interferome¬ 
ter be /i+(t, r, t, (3) and hx (t, r, f3), and let the polariza¬ 
tion angle be ^/>. Then the response of the interferometer 
will be h{t) + n(t), where n{t) is the interferometer noise, 
and 


Hi) = F+id, PKit, A L P) 

+ Fx{0,ip,'if)hx{t,r,L,P). 


(2.35) 


Here E^ and Ex are the interferometer beam pattern 
functions, given in, e.g., Ref. 0|. The dependence of the 
Fourier transformed waveform /i_|_ on r is clearly of the 
form 


h+if,r,L,P) = H+{f,L,(3)/r 


(2.36) 


for some function H+\ we define E[x {f,i,P) similarly. 
Combining Eqs. (|3), d^) and (lUI) gives 


\F+H+ + FxHx 

Shif) 


-df. (2.37) 


4 

p^{r,e,ip,il;,L,P) = — 

’’’ Jo 


We now average over all of the angles 9, ip, tf, l and (3. 
The average over polarizations and over the sky location 
gives {Fj) = {FD = 1/5, (E+Ex) = 0 0, where the 
meaning of the angular brackets is given by, for example, 

{fD = ^J dne,^£ ^F+{9, ip, (2.38) 


From Eq. (2.37) this gives 


{P^) = 


5j,2 


Hifl 

Shif) 


df, 


(2.39) 


where 


E(/)2 = ^J dn,,p + \Hx{i,P)\^) . (2.40) 

We now express the energy spectrum dE/df of the 
waves in terms of the quantity H{f)‘^. The local energy 
flux is 


dE 


1 


dA dt IGtt 


dh+ \ /dhx 


dt 


dt 


(2.41) 


where the overbar means an average over several cycles 
of the wave. Switching to the frequency domain using 
Parseval’s theorem, inserting a factor of two to account 
for the folding of negative frequencies into positive, and 
using \h+^xif)PdA = \H+^xif)pdn gives 


dE 

dVLdf 


^(|E+(6,/3)|^ + |Ex(l/3)|"). (2.42) 


Combining Eqs. ( ^.3S| ), (2.40) and ( 2.42| ) now yields 


(P^> = 


_ 

STT^r^ Jo ,/ pShif) dfldf 


df dVt- 


[/]• (2.43) 


This is essentially Eq. (2.34) with z = 0 and D{z) = r, 
the limiting form that applies when cosmological effects 
are neglected. 

Consider now so urces at cosmological distances. We 
can generalize Eq. ( ^.43 ) to inc orpo rate redshift effects as 
follows: (i) Observe that Eq. (2.43) is valid for arbitrary 
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bursts of gravitational waves provided that we interpret 
the quantity 

1 dE 


d^df 


as the locally measured energy flux dE/dAdf. (ii) Use 
the fact that the number of gravitons per unit solid an¬ 
gle per unit frequency is conserved for propagation in a 
Friedmann-Robertson-Walker background in the geomet¬ 
ric optics limit: 


dE 

d^df 


(/) = 


dEe 

dfldff 


■[(1 + ^)/]- 


(2.44) 


Here fe is the frequency at the source and / = fe /(1 + 2 ) 
is the frequency at the detector, (iii) The conversion 
factor at the detector from energy per unit solid angle to 
energy per unit area is just (1 -|- z)^/D(z)^, where D{z) 
is the luminosity distance (^. Hence 


dE 


(1 


dE^ 


■[(1 + ^)/]- 


dAdf D{zY diddfe ^ ^ 

Combining this with Eq. ( p.43 ) now yields Eq. (2.34). 

An alternati ve de rivation of the cosmological modifi¬ 
cations to Eq. ( 2.43 ), which yields a useful rule of thumb 
for understanding redshift effects, is the following. Sup¬ 
pose that an arbitrary source of waves is specified by 
parameters Oi for 1 < j < n. In the local wave zone the 
gravitational waveform is of the form 


h{t) = 




(2.46) 


for some function El, where r is the distance to the source. 
After propagating through a FRW background the wave¬ 
form becomes 


h{t) = 




redshifted 


Diz) 


(2.47) 


where if 9i has dimension (length)^’ o r (m ass)^, then 
redshifted is just (1 -|- z^Oi. Equation (2.47) can be de¬ 


rived from the general solution for wave propagation in 
FRW backgrounds in the geometric optics limit given in 
Ref. 1^. The fact that the waveform must depend on 
and reveal the “redshifted” source parameters is clear: all 
source parameters with dimension that we can measure 
from the waves must essentially come from timescales in 
the waveform, since the wave amplitude is dimensionless. 
All these timescales will be redshifted. 

There is a corresponding simple rule for how the SNR 
is affected by cosmological effects. Suppose that when 
cosmological effects are neglected, the SNR is of the form 


P = 


m) 


(2.48) 


for some function F. Then the correct SNR is just 


P = 


F{e^ 

,redshifted) 

~W) ■ 


(2.49) 


Applying this rule to Eq. ( ^43 ) yields Eq. ( ^34 ). 


III. THE GRAVITATIONAL-WAVE SIGNAL 
FROM GOALESCING BLACK HOLES 

In this section we describe our assumptions concern¬ 
ing the gravitational-wave signal from BBH mergers, and 
the evidence t hat u nderlies those assumptions. We start 
in subsection HI A by discussing the approximate split¬ 
ting of the total wavefor m into the inspiral, merger, and 
ringdown epochs. In Sec. HIB we discuss the various cir¬ 
cumstantial pieces of evidence, culled from the literature, 
about the energy spectrum and the total energy radiated 
during the nonlinear merger phase. Our assumed merger 
energy spectrum is give n by E q. ( 3.IT), with parame¬ 
ters given by Eqs. (^), ( ti.5[ ) and ( 13.131). The inspi¬ 
ral energy s pectr um is disc ussed in Sec. HI C and spec¬ 
ified in Eq. ( 3.20| ). Section HID discusses the spectrum 
and likely overall amplitude of t he ri ngdown signal; the 
spectral shape is given in Eq. ( 3.24 ) and our assump¬ 
tion abou t the overall energy radiated in Eq. ( 3.28 ). In 
Sec. HIE we estimate the SNR obtainable in band-pass 


filtering searches for waves from the merger phase. 


A. The three phases of the gravitational-wave signal 

As discussed in the Introduction, the coalescence and 
its associated gravitational-wave signal can be divided 
into three successive epochs in the time domain: inspiral, 
merger, and ringdown. Physically, the inspiral consists 
of the portion of the coalescence in which the black holes 
are separated bodies that gradually lose energy and an¬ 
gular momentum, slowly spiraling towards one another. 
By contrast, the merger consists of that portion of the 
coalescence in which the dynamics are highly nonlinear 
and must be treated by numerical relativity. With this 
in mind, it is useful to define the end of inspiral as the 
time and frequency | |7^ at which numerically generated 
templates start to be needed. Up to this time, post- 
Newtonian templates, possibly s upp lemented with IBBH 
templates, will be used (c/. Sec. IB). 

After inspiral, we expect there to be a lot of compli¬ 
cated merger dynamics that at present are not well un¬ 
derstood. The system will gradually settle down to a 
Kerr black hole; the last gravitational waves we expect 
to see are those produced by the quasi-normal ringing 
modes of this merged black hole. It is clear that there 
will be a smooth transition in the gravitational wave¬ 
form from the merger portion to the ringdown portion, 
as the effects of nonlinearities become less and less im¬ 
portant with time. As this happens, the signal should 
become increasingly well approximated by a linear com¬ 
bination of exponentially decaying sine waves. This is 
the behavior that has been seen in numerical simulations 
of, for example, head-on collisions ||7^ , [^ . At late times, 
the I = m = 2 mode will dominate over other quasi¬ 
normal modes, for two reasons which are of comparable 
importance: (i) The I = m = 2 mode is the most slowly 
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(3.1) 


damped of all the QNR modes Q, and (ii) during co¬ 
alescence, the binary will have a rotating shape roughly 
corresponding to spheroidal harmonic indices I = m = 2, 
and thus this mode will be preferentially excited [^. 
We define the ringdown as beginning when the wave¬ 
form becomes dominated by the I = m = 2 QNR mode; 
the merger thus contains those portions of the waveform 
where other modes and/or non-linear mode-mode cou¬ 
plings are important. Clearly there is some arbitrariness 
in the exact time at which the ringdown starts, related 
to the accuracy we require of the fit of the waveform to 
the ringdown signal. 

Note that, by definition, the three phases of the sig¬ 
nal are disjoint in the time domain. It does not follow, 
though, that they should be disjoint in the frequency do¬ 
main: their energy spectra might overlap. However, it is 
at least approximately true that the inspiral and merger 
and disjoint in both time and frequency. Let tmerge be the 
time at which inspiral is said to end and merger to be¬ 
gin. Then, since the inspiral waves chirp monotonically 
and since the adiabatic approximation is just beginning 
to break down at the end of the inspiral, the inspiral 
spectrum is confined, to good approximation, to the fre¬ 
quency region / < /merge- Here, /merge is the frequency 
of the inspiral waves at t = tmerge- We shall further as¬ 
sume that the merger waves’ spectrum is confined to the 
frequency regime / > /merge; this assumption should be 
valid to a moderately good approximation. (We discuss 
below estimates of the value of /merge-) 

One particular component of the gravitational-wave 
signal, the Christodoulou memory [7^ , will violate 
this assumption. The memory component of the 
gravitational-wave signal has most of its power below 
/merge in the frequency domain, but accumulates gradu¬ 
ally during the inspiral, merger and ringdown in the time 
domain. The memory waves from BBH coalescences will 
probably not be detectable with ground based interfer¬ 
ometers, but very probably will be detectable with LISA 
[fz^ . We will neglect the memory component of the waves 
in our analysis, since it will not be as easy to detect as 
the components we do discuss. 


B. Energy spectrum of the emitted gravitational 
radiation from the merger phase 

The total amount of energy radiated in BBH merg¬ 
ers, and its distribution in frequency, is highly uncertain 
because detailed numerical calculations of these mergers 
have not yet been made. In this subsection, we discuss 
what little evidence there is about the energy radiated, 
and describe our crude model of the spectrum. We as¬ 
sume a uniform distribution in frequency from the lower 
frequency /merge mentioned above and some upper fre¬ 
quency /high- 

The total amount of energy radiated during a BBH 
coalescence will be some fraction e of the total mass M = 
nil + m 2 of the system: 


Aradiated — cAf. 


The fraction e will depend only on the mass ratio mi/m 2 , 
on the initial spins Si and S 2 of the two black holes, and 
on the initial direction L of the orbital angular momen¬ 
tum ra: 


mi 


r,L 


M2’ M2’ 

We can very roughly divide up this fraction as 

e — Cinspiral T ^merger T Cringdown, 


(3.2) 


(3.3) 


according to the amounts of energy radiated in the three 
different epochs of the waveform. There is clearly some 
arbitrariness in this division of the radiated energy into 
three pieces, which is related to the choice of the fre¬ 
quency /merge at which the inspiral is said to end and 
the merger to begin |Q, and also to the choice of the 
time tqnr at which the merger is said to end and the ring- 
down to begin. 

We now discuss methods of estimating /merge, Cmerger, 
and Cringdo-wn- Roughly speaking, we want to choose the 
interface between “inspiral” and “merger” to be at that 
point where post-Newtonian templates cease to be useful 
and where numerically generated templates will need to 
be used. Numerically generated initial data sets of two 
orbiting black holes provide some guidance about this 
interface. Cook compares the properties of such initial 
data sets to the predictions of post-Newtonian theory 
at second post-Newtonian order. He finds, for example, 
that the discrepancy in the binding energy will be ~ 5% 
when the gravitational-wave frequency (twice the orbital 
frequency) is / ^ 0.02/M, where M is the total mass of 
the system, and will be ^ 15% at / ^ 0.05/A/ ||^. 

Towards the end of the inspiral, the black holes make a 
gradual transition from a radiation-reaction driven inspi¬ 
ral to a freely-falling plunge after which, even 

if the radiation reaction could be turned off, the black 
holes would still merge. This will occur roughly at the 
last stable circular orbit (LSCO). It would seem to make 
sense to choose /merge = /lsco, the frequency at which 
this occurs. Strictly speaking, the concept of an LSCO 
makes sense only in the test particle limit (mi/m 2 <SC I). 
Various somewhat ad hoc methods have been used to es¬ 
timate /lsco outside this regime, but the results differ 
by factors ^ 2. Cook’s initial data set analysis together 
with the calculation of an “effective potential” yields the 
estimate /lsco 0.055/M for equal mass black holes 
1^. In post-Newtonian theory, the LSCO can be defined 
by artificially turning off the radiation reaction terms in 
the equations of motion. Using this method, Kidder, 
Will and Wiseman estimate that /lsco ~ 0.02/M [|j. 
They use hybrid equations of motion which are accurate 
to post-2-Newtonian order for comparable mass ratios, 
and which are exact in the test particle limit. The value 
they obtain varies by less than ^ 20% as the mass ra¬ 
tio is varied. Note that the frequency 0.02/M is close 
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to the value 6“^/^/(7rM) obtained from the usual LSCO 
at r = 6M in the test particle limit (here r is the usual 
Schwarzschild radial coordinate). Finally, earlier analy¬ 
ses by Blackburn and Detweiler used a variational prin¬ 
ciple together with the assumption of periodic solutions 
to Einstein’s equations to obtain the approximate lower 
bound /lsco ^ 0.06/M |^|. 

Note that all of these estimates are for equal mass, non 
spinning black holes; the value of the frequency /lsco 
can also vary by factors ^ 2 if the black holes are spinning 
or have different masses. 

In light of this uncertainty, we somewhat arbitrarily 
adopt the conservative value of 


/. 


merge 


0.02 

'W 


= 205 Hz 


/2OM0 
I M 


(3.4) 


as the lower cutoff frequency for the energy spectrum. 
This (low) value of /merge is conservative in the sense 
that we can be reasonably sure that numerically gener¬ 
ated templates will not be needed before / = /merge- On 
the other hand, it may optimistically overestimate the 
merger SNR by increasing the number of cycles in what 
we define as our merger waveform at the expense of the 
number of cycles in the inspiral. 

We next discuss our choice of upper frequency shut¬ 
off for the merger energy spectrum. As discussed above, 
we define the end of the merger to occur at that time 
tqnr after which the waveform can be fit fairly accurately 
by the I = m = 2 QNR mode signal. The merger and 
the ringdown will therefore be disjoint in the time do¬ 
main, but not necessarily in the frequency domain. It 
seems likely, however, that an approximate upper bound 
for the frequencies carrying appreciable power during the 
merger is the quasinormal ringing frequency itself. This 
is supported to some extent by calculations of the energy 
spectrum of waves produced in the test particle limit— 
see Fig. 2 of Ref. |^|. The energy spectrum dE/df in 
that case is peaked near the quasinormal mode frequency 
.06/M of a non-rotating black hole, and has fallen by 2 
orders of magnitude at .14/M. Although this analysis 
does not incorporate any nonlinear effects, there is some 
evidence that the energy spectrum in the equal mass 
case will be qualitatively similar—namely, the head-on 
collision of black holes in full numerical relativity has 
been computed, and the resulting waveforms are qualita- 
tively similar to those predicted from perturbation the¬ 
ory It is not clear how relevant the head-on col¬ 

lision case is to the gradual inspiral; but on the other 
hand, there is no other guidance available at the present 
time. We will therefore use as our upper shutoff fre¬ 
quency for the merger spectrum the frequency /qnr of 
the I = m = 2 quasinormal mode. This frequency will 
depend on the unknown dimensionless spin parameter a 
of the final Kerr black hole (its spin angular momentum 
divided by M^). We choose to use the value of a for 


which the frequency is highest: in the limit a —> I, the 
QNR frequency becomes 0.9/(27rM) (more than twice 
the value for a Schwarzschild black hole) |Q. Thus, our 
upper cutoff frequency is 


/high — /qnr — 


0.13 


= 1430 Hz 


/ 20Me \ 

\ M )■ 


(3.5) 


Our re asons for assuming a high value of a are discussed 
in Sec. HID below. 

Finally, consider the total amount of energy 


Aj-ad (^merger T ^ringdown )-^1^ 


(3.6) 


radiated during the final merger and ringdown. To esti¬ 
mate this we consider first an extrapolation due to Smarr 
@0- The idea is to calculate the energy radiated in the 
test particle approximation using black-hole perturbation 
theory. The result will be of the form Arad = fc/r^/M, 
where A: is a dimensionless constant, /r is the mass of the 
test particle, and M is the total mass. Now simply re¬ 
place ^ by the reduced mass mim 2 /(mi -I-TO 2 ) of the two 
black holes. This extrapolation works exactly in Newto¬ 
nian gravity supplemented by the quadrupole formula, in 
the sense that the dynamics of two bodies of comparable 
masses can be deduced from the corresponding dynam¬ 
ics in the extreme mass ratio limit. Surprisingly, it also 
works to within ~ 20% to predict the total energy ra¬ 
diated in the head-on collision of two equal-mass black 
holes Thus, one might consider applying 

this extrapolation to inspiral-preceded mergers. Such a 
calculation would require a determination of the radiated 
energy from the final plunge of a test particle, he., the 
solution of a coupled radiation reaction/radiation gener¬ 
ation problem outside of the adiabatic regime (although 
the coupling is important only in a brief transitional 
phase and essentially sets the initial conditions for an 
approximately geodesic plunge [^). This calculation 
has not yet been done. However, there is an alterna¬ 
tive approximate method to estimate the amount of en¬ 
ergy radiated in the test particle limit. The innermost, 
marginally bound (unstable) circular orbits should be a 
fair approximation to the motion during the final plunge. 
Detweiler 0| showed that the amount of energy radiated 
per orbit by a test particle in such an orbit is of the form 
Arad = kfjL^/M, where k varies from 0.65 at a = 0 to 2.8 
at a = 0.95, where a is the spin parameter of the black 
hole. Assuming that there will be ^ 1 effective orbit 
during the final plunge, Detweiler estimated that [ 0 

0mMF{fi/M) < Arad ^ 0.2M A(^/M), (3.7) 

where 

Fid/M) = (3.8) 

which is unity in the equal-mass case /r = M/4. 
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Arguments based on conservation of angular momen¬ 
tum also suggest an approximate lower bound on the 
radiated energy of about O.IM for equal-mass BBHs in 
the most favorable cases, as we now outline. Roughly 
speaking, the system’s angular momentum divides up as 

Si -f S 2 -|- Lorb = Lrad + Sfinal, (3-9) 

where Si and S 2 are the spins of the two black holes 
just before the final plunge, Lorb is the orbital angular 
momentum just before the plunge, Lrad is the angular 
momentum carried off by the gravitational waves, and 
Sfinal is the spin of the final Kerr black hole. This split¬ 
ting of the total angular momentum of the spacetime into 
various pieces is really only well defined only in a post- 
Newtonian type of limit; however, the effects of this am¬ 
biguity are presumably not important for the purposes of 
our crude, approximate argument. We now specialize to 
the most favorable case in which Si, S 2 and Lorb are all 
aligned. We take the orbital angular momentum at the 
beginning of merger to be |Lorb| ~ 0.9M^, the value pre¬ 
dicted by Cook’s initial data sets at / = 0.02/M ||]. We 
also assume that both inspiraling black holes are rapidly 
spinning, so that |Si| « IS 2 I « (M/2)^. Equation ( |3.9| ) 
then yields 

ILradI >0.4M2, (3.10) 

since |Sfi„ai| < M^. 

Next, we use that fact that the energy E^ad and the 
angular momentum Lrad carried off by gravitons of fre¬ 
quency / and azimuthal multipole order m are related 

by [H 


Lirad — f Tfi. 


(3.11) 


If we estimate / as (/merge + /qnr)/2, and make the admit¬ 
tedly optimistic assumption that m ost of the ra diati on is 
quadrupolar, we obtain from Eqs. ( 3.10 ) and ( 3.11 ) the 
crude estimate for the total radiated energy M 


Lrad > O.IM. 


(3.12) 


This estimate includes both merger and ringdown radi¬ 
ated energies; we need to subtract out the ringdown por¬ 
tion to obtain the energy contained in the merger portion 
of the waveform. Below we estimate ~ 0.03M to be an 
approximate upper boun d for the ringdown energy, and 
hence most of the energy ( 3.12 ) should be radiated in the 
merger waves. 

There is an additional, separate argument on e can 
make which indicates that most of the energy ( 3.12| ) 
should be radiated as merger waves and not as ring- 
down waves. As noted by Eardley and Hirschmann [p4[ , 
any system with J > M^ cannot evolve to J < M^ 
by radiating quadrupolar waves at the ringing frequency 
/qnr ~ 1/{2'kM) of a near-cxtremal Kerr black hole. This 
is because at this high frequency, too much mass-ene rgy i s 
radiated per unit angular momentum radiated; Eq. (3.11) 


with m = 2 and with / = /q^r yields AJ = A(M^). 
Hence, since the final black hole must have J < M^, a 
substantial amount of the radiation must be emitted at 
lower frequencies. 


Based on the estimates (3.7) and ( 3.12 ), and on the 
estimated upper bound ~ 0.03M mentioned above for 
the ringdown radiated energy, we perhaps optimistically 
take O.IM to be the energy radiated during merger in the 
equal-mass case (which corresponds to about 0.13M total 
radiated energy in the merger and ringdown). For non 
equal-mass BBHs we a ssum e that the radiated energy is 
reduced by the factor (|3.8|), so that 


L 


merger 


^merger 


F{fi/M)M = 0.1F{fi/M)M. (3.13) 


The rather high radiation efficiency of 0.1 that we are 
assuming is probably most plausible in the context of 
rapidly spinning coalescing black holes. In particular, if 
the spins and the orbital angular momentum are some¬ 
what misaligned, the coalescence will be something like 
“two tornados with orientations skewed to each other, 
embedded in a third, larger tornado with a third ori¬ 
entation” j^. Intuitively one would expect that such 
systems have more “settling down” to do to get to the 
final Kerr black hole, and that correspondingly the non¬ 
linear, highly dynamical phase should last longer and/or 
produce more radiation. Also, the potential barrier that 
surrounds the final black hole (which normally tends to 
reflect back into the black hole the dominant waves of 
frequency / ~ 1/a few times M) presumably will effec¬ 
tively not be present during the violent phase of a merger 
in which the spins and orbital angular momentum are of 
comparable magnitude and are misaligned. 

We note that coalescences which radiate as much en¬ 


ergy as the estimate (3.13) will also likely radiate a sub¬ 
stantial amount of linear momentum, and the consequent 
recoil of the final black hole could correspond to a kick 
velocity that is a moderate fraction of the speed of light. 

Finally, consider the shape of the energy spectrum 
dE/df between /merge and /qnr- For simplicity, and for 
lack of evidence in favor of anything more specific, we 
choose a flat spectrum: 


dE €mMF{fi/M) 

df /qnr /merge 

= 0mM^F{^l/M) 0 (/ - /. 


©(/ - /merge)0(/qnr “ /) 

merge)0(/qnr " /), (3.14) 


where Cm = Cmerger = 0.1 and 0 is the step function. 

Our assumed merger spectrum has the desirable fea¬ 
ture that it is reasonably consistent with the inspiral en¬ 
ergy spectrum at / = /merge- This can be seen as follows. 
We do not expect a very large change in dE/df at the 
interface between inspiral and merger, for the following 
reason: throughout the inspiral, the energy spectrum can 
be written 

dR 

— (/)cxn,,,(/)/^[^(/)]^ (3.15) 
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where ncydf) is the number of cycles spent near fre¬ 
quency /, and h[t{f)] is the wave amplitude at the time 
t{f) when the waves’ frequency is / [see also Eqs. (2.£) 
and (5.2)]. In the transition from inspiral to plunge, 
the amplitude h[t{f)] should change only moderately, 
but ricycif) could conceivably change by a large amount: 
from its pre-plunge value to ^ 1. It turns out that the 
number of cycles does not change very drastically: at the 
interface between inspiral and merger. 


= (/) = 


Z! 

/ 

6.7 


9071 


(irMf) 


1 -5/3 


-5/3 


/ merge 


(3.16) 


(Here, M = is the so-called chirp mass.) 

Hence, we do not expect dE/df to change too much. 

We now show that our assumed values of dE/df from 
the end of the inspiral and from the beginning of the 
merger are reasonably close to one another, in agreement 
with the above expectation. Using the approximation 
that the motion is an adiabatic evolution from one circu¬ 
lar orbit to another, the energy radiated during inspiral 
can be written 


w) 


radiated 


= ^^^circ(/) [1 + 0(l/neyc(/))] , (3.17) 


where Edrdf) is the binding energy of a particle in 
a circular orbit at orbital frequency f/2 (recall that 
the gravitational-wave frequency is twice the orbital fre¬ 
quency). At our estimated plunge frequency /merge = 
.02/M, the hybrid equations of motion of Kidder, Will 
and Wiseman predict in the equal mass case that 


dEc 


-(/ = 0.02/M) = 0.88M^ 


df 

while Cook’s initial data sets predict that 
dE ■ 

^-^circ , 


df 


-(/ = 0.02/M) = 0.41MC 


(3.18) 


(3.19) 


The true value of dEdrc/df is presumably somewhere be¬ 
tween these two values. We would expect the value of 
((i£’/(i/)radiated_to be Smaller than this, because the ap¬ 
proximation (3.17) is breaking down at this point—some 
of the gravitational binding energy is being converted 
to inward radial kinetic energy instead of being radiated. 
This picture is roughly consistent with our assumpt ion o f 
dE/ df = 0 .91M^ for / > 0.02/M, given that Eqs. (ItTI) 
and ( |3.19| ) apply to the non-spinning case, and we are 
interested in the rapidly spinning case for which dE/df 
could be a few times larger. 


C. Energy spectrum of the radiation from the 
inspiral phase 

The standard quadrupole formula prediction for the 
inspiral energy spectrum is (see, e.g., Ref. [^) 


^ = i^2/3 M2/3^-1/3^ 

df 3 


(3.20) 


This approximate formula is adequate to estimate the 
SNR obtained from optimal filtering of the inspiral wave¬ 
form. The formula will be accurate to within a few 
tens of percent up to / = /merge @. Using Eq. ( 3.20| ) 
to estimate the SNR effectively assumes that both the 
gravitational-wave signal and the theoretical templates 
used to cross-correlate with the d ata st ream are given 
by the quadrupole approximation ( ^.20 ). The SNR we 
calculate using this lowest order formula will be approx¬ 
imately the same as that found by cross-correlating real 
sign als [w hich incorporate higher order corrections to 
Eq. ( 3.20| )] against sufficiently accurate theoretical tem¬ 
plates. As outlined in Sec. [B, the required template 
accuracy should be achievable by post-Newtonian ex¬ 
pansions 1^1^, perhaps supplemented with alternative 


techniques for the latter, high frequency part of the sig¬ 
nal at 0.01/M d f d 0.02/M (the IBBH regime). We 
assume that the inspiral energy spectru m shu ts off at 
/ = /merge = 0.02/M, as discusscd in Sec. HI B|. 


D. Energy spectrum of the radiation from the 
ringdown phase 


As we have discussed, the ringdown portion of the 
gravitational-wave signal is that portion in the time do¬ 
main that can be fit fairly accurately by an exponen¬ 
tially decaying sinusoid corresponding to the most slowly 
damped I = m = 2 quasinormal mode of the final black 
hole. The shape of the corresponding energy spectrum is 
well understood: it is just a resonance curve (although 
see Appendix B for discussion of a subtlety in the ap¬ 
plicability of the concept of the waves’ energy spectrum 
for calculating ringdown SNRs). The overall amplitude 
of the energy spectrum, however, is not well understood. 
In this section we discuss our estimate of the overall am¬ 
plitude for the ringdown signal. 

The QNR gravitational waveforms h+ ft, l, /3) and 
/lx it, L, (3) are given by @ 


/i_|_ — ihx 


-2<52(bP,a) 

r 

X exp[-2/7r/qnrt “ I/t + 


(3.21) 


for t > 0. Here we have chosen t = 0 to be the time of the 
start of the ringdown, M is the mass of the black hole, 
oM^ is its spin, and (po is a constant phase. The quanti¬ 
ties t and jS are sphe rical polar coordinates centered on 
the black hole [cf. Sec. H C |, 2 'S'|(t, /3, a) is a spin weighted 
spheroidal harmonic whose angle averaged rms value is 



(3.22) 


and xl is a dimensionless coefficient that describes the 
magnitude of the perturbation when the ringdown be¬ 
gins. The quantities /qnr and r are the frequency and 
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damping time, respectively, of the QNR mode of the 
black hole. The quality factor Q of the mode is given 

by Q = TTT/qnr- 

As mentioned in the Introduction, there is a mapping 
between the values of M and a, and of fnnr and r. This 
mappiM has been explored by Leaver |g7|] and Echev- 
erria Leaver studied in depth the damping times 

and ringing frequencies of various QNR modes for var¬ 
ious Kerr black hole masses and spins. By numerically 
solving the Teukolsky equation, he produced extensive 
catalogs of the variation of r and with black hole 

mass M and dimensionless spin a From that data 

for the I = m = 2 mode, Echeverria produced the 
following analytic fits, which are good to within roughly 
5%: 


/qnr 


Q 


1-0.63(1 
1-0.63(1 
2(l-a)“®/^°. 


1 


2ttM 
f 37MQ \ 
\ M ) 


865 Hz 

(3.23) 


Note that Q diverges as a ^ 1 but that /qnr tends to¬ 
wards a constant. 

The energy spectrum for the QNR waveform ( 3.2l| ) is 
derived in Appendix B and is given by 


dE _ f 1 

If ~ 327r3 r2 \ [(/ _ + (2^^)-2]2 

+- - - 2 I (3-24) 

[(/ + /qnr)2 + (2^T)-2f / 

« i^'QM'/qrrr <5(/ " /qnr) [1 + 0(1/Q)] ■ (3.25) 


mass BBHs in the following way. When the I = m = 2 
quasinormal mode becomes the dominant mode present, 
the merged body will be a distorted Kerr black hole. The 
horizon’s cross section, seen looking down on this black 
hole from above, will be a rotating oval rather than a 
circle. We can quantify the distortion by the ratio of the 
polar circumference about the long axis of this oval to 
the polar circumference about the short axis. Let A 2 de¬ 
note the perturbation amplitude such that the ratio of 
circumferences is, say, 2:1. It is fairly clear that non¬ 
linear couplings between the I = m = 2 mode and other 
modes will be important for A -f, A 2 ; thus, at this 2 : 1 
distortion ratio , the signal will not be well approximated 
by just the I = m = 2 mode. Hence, A 2 is a reasonable 
upper bound for the true amplitude A. 

We could in principle calculate the upper bound A 2 by 
writing the spacetime metric as 


„ _ „KERR 

9ab — 9ab 


A 2 h 


QNR 

ab 


where is the Kerr metric and I = m = 

2 quasinormal mod e whose asymptotic form at large r is 
given by Eq. ( ^.21 ), and by calculating from this met¬ 
ric the ratio of circumferences [^. For this paper, we 
used a much less sophisticated method to estimate A 2 . 
Using the quadrupole formula, we examined the radia¬ 
tion produced by a solid body that is distorted with the 
same 2 : 1 circumference ratio as the merged black hole, 
and obtained the estimate A 2 ~ 0.4 |Q. Setting our 
waveform amplitude A equal to this upper bound yields 
an rms angle-averaged waveform h = (0.4/V^)( M/r) = 
0.1(M/r) at the beginning of ringdown. From Eq. ( ^.24 ), 
the corresponding radiated energy is 


Aringdown ~ « 0.03M. 


(3.26) 


Approxima ting the energy spectrum by a delta function 
as in Eq. ( 3.25 ) will often (but not always) provide a 
fairly good approximation to the SNR; see Appendix B 
for more details. 

The values of the spin a of the final black hole (and 
correspondingly, Q) and also of the amplitude A will 
depend on the initial parameters of the system as in 
Eq. (3.2). This dependence is very poorly understood 
at present. We expect the final bl ack ho le to be rapidly 
spinning since, as explained in Sec. HIB, the total angu¬ 
lar momentum of the binary at the end of the inspiral is 
^ 0.9M^ when even the individual black holes are non¬ 
spinning 1^^ , and the individual black hole spins can add 
to this. Moreover, the individual black holes in the bi¬ 
nary may typically have been spun up to near maximal 
rotation by an accretion disk p8|] . For definiteness, we 
somewhat arbi trarily take a = 0.98, which corresponds 
from Eq. ( 3.23|) to Q = 12 and /qnr = 0.13/M. The final 
ringdown SNR values we obtain vary only weakly with 
our assumed value of a [cf Eq. ( A14 )], for fixed total 
energy radiated in the ringdown. 

Although the value of the overall amplitude A is un¬ 
certain, we can estimate a upper bound on A for equal 


As mentioned in the Introduction, comparable ringdown 
radiation efficiencies ~ 3% have been seen in numerical 
simulations of the evolution of distorted, spinning black 
holes iQ. 

To summarize, our assumed parameter values for the 
black hole dimensionless spin parameter a and for the di¬ 
mensionless amplitude parameter A for equal-mass BBHs 
are 


a = 0.98 

A = 0.4. (3.27) 


These imply the following values for the quasi-normal 
ringing frequency /qnr, the quality factor Q and the 
fraction of total mass-energy radiated during ringdown 

Cringdown — jngdown / Af. 


_ 0.13 

/qnr- ^ 
12 


^ringdown 


0.03. 


(3.28) 


For non equal-mass BBHs, we a ssu me that A and 
Cringdown ars reduced by the factor (3.8). 
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E. Number of independent frequency bins for the 
merger phase 


In Sec. IIA we showed that for any burst of gravita¬ 


tional waves, the band-pass filtering SNR is smaller than 
the optimal filtering SNR by a factor of approximately 


= V27’A/. 


(3.29) 


Here T is the duration of the merger and A/ is the 
band-pass filter’s bandwidth, which should be approxi¬ 
mately the largest expected signal bandwidth. The prod¬ 
uct JVhins, which is roughly the number of independent 
frequency bins (or Fourier coefficients) needed to describe 
the signal, also determines how well matched filtering 
searches perform com pare d to noise-monitoring searches, 
as explained in Sec. OH. In this section, we obtain a 
rough estimate of the parameter A/bins for the merger 
gravitational waves. 

First consider the bandwidth A/. Our assumed band¬ 
width for the merger signal is A/ = /qnr - /merge ~ /qnr 
[since /merge < /qnr; c/. Eqs. (M ) and dU)]. We cannot, 
however, we completely confident that all signal power 
in the merger will lie at frequencies below /qnr, so a 
more appropriate choice might be A/ ~ 2/q„r. Also 
the quasinormal ringing frequency /qnr depends on the 
dimensionless spin parameter a of the final black hole as 
given by Eq. ( |3.23| ). Choosing the highest possible value, 
/qnr = l/(27rM)], yields 

(3-30) 


Turn, now, to the effective duration T of the signal, 
which is defined by Eq. ( 2.13| ). Clearly this parameter 
will vary considerably from event to event and its value is 
highly uncertain. To get a feeling for the range of possible 
valu es, let us consider the type of coalescence discussed in 
Sec. IIIB , where both inspiraling black holes are nearly 
maximally spinning, with their spins and the orbital an¬ 
gular momentum of the bina ry all aligned. In this favor¬ 
able case, as argued in Sec. IIIB, at the end of inspiral 
near r ~ 6M the binary has an excess ~ 0.4M^ of total 
angular momentum that it needs to shed before it can 
settle down to its final Kerr state. Thus it seems reason¬ 
able to suppose that the two black holes will be somewhat 
centrifugally hung up and continue to orbit for several cy¬ 
cles before their event horizons merge. In other words, 
the orbital dynamical instability may be suppressed by 
the aligned spins, so that what we call merger (for data- 
analysis and computational purposes) may dynamically 
resemble the inspiral. The duration of the merger in this 
case may therefore be long. By contrast, for inspiral¬ 
ing non-spinning black holes, there is probably no excess 
angular momentum that needs to be radiated after the 
orbital dynamical instability is reached, and so it is plau¬ 
sible that the two black holes merge rather quickly. In 
such a case, most of the emitted energy could come out 


in the ringdown waves rather than in merger waves, as 
has typically been seen in numerical simulations. 

To estimate the duration T for the first type of sce¬ 
nario, we assume that the gravitational-wave luminosity 
dE/dt during the merger phase is roughly the same as the 
luminosity at the start of the ringdown, 2erM/T. (Here, 
er = 0.03 is the ringdown radiation efficiency, and r is 
the damping time of the QNR mode). Combining this 
with a total energy radiated during the merger of emM 
(with Cm = 0.1) yields the estimate 

T=i^r. (3.31) 


Clearly this estimate will become invalid for high val¬ 
ues of r; in that limit, the high quality factor of the 
QNR mode causes a low quasi-normal ringing luminos¬ 
ity, whereas there is no reason for the merger luminosity 
to be comparably low. Nevertheless, we insert our as¬ 
sumed value of r , wh ich cor respo nds to a q uality factor 
of Q = 12 [Eqs. ( 3.23 ) and ( [1.27 )] into Eq. (3.31) to ob¬ 
tain T ~ 50M. Combining this with Eqs. ( ^.14 ) and 
yields the estimate 


v'AAbins - 730 - 5. 


(3.32) 


Eor inspiraling Schwarzschild black holes, on the other 
hand, T may not be much larger than r = lOM (assum¬ 
ing a = 0.5 say) , givin g VA^bins ^ VG- 

The factor VA^bins is thus likely to lie in the range 
2 ^ y/Ab ins ^ 5. We adopt the estimate VAbins = 4 in 
Sec. VIB below to estimate the reduction in SNR result¬ 
ing from using band-pass filters instead of optimal tem¬ 
plates. We u se the conservative, large value Abins = 60 
in Sec. VIB to estimate detection thresholds for noise¬ 
monitoring searches for signals. 


IV. INTERFEROMETER NOISE CURVES 


In this section we describe our approximate, piecewise 
power law, analytic model for the noise curves for the 
initial LIGO interferometers, for the advanced LIGO in¬ 
terferometers, and for the LISA interferometer. We ex¬ 
press our model in terms of the dimensionless quantity 
hrmsif) = VfShif), where Sh{f) is the one sided power 
spectral density of the interferometer noise |Q. The 
quantity /irms(/) is the rms fluctuation in the noise at 
frequency / in a bandwidth A/ = / 0. Our approxi¬ 
mate model for the noise spectrum is 

f oo f < fs, 

h (f)- } (a///m)”^/^ fs< f < fm/a 

- S U/a <f<afm 

[ hm afm < /■ 

This noise curve scales like /^A at high frequencies and 
like /“^A at low frequencies, and has a flat portion at in¬ 
termediate frequencies. The noise curve depends on four 
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parameters: (i) A lower shutoff frequency fs below which 
the interferometer noise rapidly becomes very large and 
can be approximated to be infinite. For the ground based 
interferometers, this low-frequency shutoff is due to seis¬ 
mic noise, while for LISA it is due to accelerometer noise 
(Ref. 1^, p. 23). (ii) A frequency fm, which is the lo¬ 
cation of the center of the flat portion of the spectrum, 
(iii) A dimensionless parameter hm, which is the mini¬ 
mum value of h^msif)- (iv) A dimensionless parameter a 
which determines the width in frequency of the flat por¬ 
tion of the noise curve. We approximate the noise curves 
by piecewise power laws in this way for calculational con¬ 
venience. 

For the initial and advanced LIGO interferometers, we 
determined best-fit values of the parameters fs, fm, hm 
and a by htting to the noise curves given in Ref. [^. 
(Note that Fig. 7 of Ref. is a factor of 3 too small 
from ~ 10 Hz to ~ 70 Hz. This error does not appear in 
Fig. 10 of that reference The resulting parameter 

values are: 

INITIAL LIGO , . 

INTERFEROMETER, 


ADVANCED LIGO , , 

INTERFEROMETER. 


fs = 40HZ ^ 

fm = 160 Hz \ 

a = lA t 

hm = 3.1 X 10-22 J 

and 

fs = lOHz 
fm = 68 Hz 
a = 1.6 

= 1.4 X 10-23 ^ 


For ground-based interfero met ers, the / 3/2 portion of 
our approximate formula Q is the thermal suspen¬ 
sion noise and the /3/2 portion models the laser shot 
noise. The intermediate portion smoothes the transition 


between these two pieces of the spectrum \M 


For the space-based LISA interferometer, we deter¬ 
mined best-fit values of the parameters fm, hm and a 
by htting to the noise curve given in Ref. [Q, and ob¬ 
tained the lower cutoff frequency fs from Ref. |^. The 
resulting parameter values are: 


u 

u 

a 

hn 


10-“^ Hz 

3.7 X 10-3 Hz 
5.5 

5.8 X 10-22 


LISA 

INTERFEROMETER. 


(4.4) 


Our piecewise power-law model Hts the true noise curve 
less accurately for LISA than for the LIGO interferome¬ 
ters, but it is still a fairly good approximation. 

The sensitivity of LISA at the lower end of its fre¬ 
quency window may be degraded somewhat by a back¬ 
ground of gravitational waves from white dwarf binaries 
0. We neglect this issue here as this white dwarf noise 
level is fairly uncertain. 


V. SIGNAL-TO-NOISE RATIOS 


In this section we calculate the angle-averaged SNRs 
for the three different stages of the coalescence (inspiral, 
merger, and ringdown) for initial LIGO interferometers, 
for advanced LIGO interferometers, and for LISA. We 
do so for black hole binaries at various distances and 
of various masses. These SNR valu es are o btain ed by 
substituting the e nerg y spectra ( 3.2C ), ( 3.14 ) and ( 3.25| ) 
estimated in Sec . |Hlt togeth er w ith the interferometer 
noise curves (4.2), ( ^!)^) an d (4^), into the SNR formula 
derived in Sec. |II Q. The analytic calculations are 
carried out in Appendix A; in this section we summarize 
the results in three graphs (Figs. ^ ||and|^, and discuss 
their implications. 

We start in subsection V A by graphing the waves’ en¬ 
ergy spectra and the interferometers’ noise spectra for 
several different examples of black hole binaries, in order 
to illustrate the factors determining the SN R va lues. The 
general results are presented in subsection VB. 


A. Specific examples 


We start by rewriting the general formula ( |2.34 ) for the 
SNR in a more useful form. If we define the characteristic 
gravitational-wave amplitude 


/ichar(/)" - [(1 




(5.1) 


then from Sec. 11C it can be seen that the SNR squared 
for an optimally oriented source is 


doptimal orientation 


J rf(ln/) 


hcharifY 

hrms(/)2 


(5.2) 


where hrms(/) = \/fShif)- However, Eq. ( 2.34 ) shows 
that the angle-averaged SNR squared is a factor of 5 
smaller than the optimal value (5.2), so we can rewrite 
Eq. (2.34) as 


(P^> = J d{lnf) 


^char(/)^ 

hnifr 


(5.3) 


where hn{f) = V^hrmsif)- We interpret the quantity 
hn as the rms noise appropriate for waves from random 
directions with random orientations |Q. The useful¬ 
ness of the quantities hchar(/) and h„(/) is that plotting 
/ichar(/) and hn{f) for various sources illustrates [from 
Eq. (^)] the possible SNR values and the distribution 
of SNR squared with frequency (see, e.g., Ref. 0). 
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FIG. 1. An illustration of the relative magnitudes of our 
estimates of the inspiral, merger and ringdown energy spec¬ 
tra in two different ca ses. The solid lines are the rms noise 
amplitudes hn{f) = \/5fSh{f) for our assumed model ( 0 ) 
of the LIGO initial and advanced interferometer noise spec¬ 
tra. The dashed and dotte d lines show the characteris¬ 
tic amplitude /ichar(/) oc ^JdE/df of the waves, defined 
by Eq. (5.1). The dehnition of /ichar is such that the sig- 
nal-to-noise ratio squared for a randomly oriented source is 
given by {S/N)'^ = J d(lnf) [hchar(/)/hn(/)]^ . The upper 
dashed and dotted lines correspond to a binary of two IOMq 
black holes, coalescing at a distance of D = 200 Mpc. The 
sloped portion of the dashed line is the inspiral part of the en¬ 
ergy spectrum, which gives an SNR for the initial (advanced) 
interferometer noise curve of 2.6 (84). The flat portion is 
our crude model of the merger part of the energy spectrum, 
which gives an SNR of 2.1 (16). The dotted line is our esti¬ 
mate of the energy spectrum from the ringdown, which gives 
an SNR of 0.1 (0.86). An event rate of very roughly 3yr“^ is 
expected for mergers inside this distance [[17|Jlq |. The lower 
dashed and dotted lines correspond to a binary of two 15Mq 
black holes at a redshift of 2 = 1 (or at a luminosity dis¬ 
tance of D = 4.6 Gpc; the cosmological parameters flo = 1 
and Ho = 75kms“^Mpc“^ were assumed). In this case the 
inspiral, merger and ringdown SNR values for the initial (ad¬ 
vanced) interferometers are 0.08, 0.42, and 0.07 (6.6, 7.2 and 
0.5) respectively. Black hole binaries with constituents this 
massive will be visible to great distances, making them a pos¬ 
sible important source, depending on the very uncertain event 
rate. The signal-to-noise ratio from the merger phase is en¬ 
hanced for these massive distant sources, in part because the 
combination of cosmological redshift plus lower intrinsic fre¬ 
quency (due to higher mass) brings the merger part of the 
energy spectrum down to lower frequencies at which the in¬ 
terferometer noise is smaller. 


In Fig. we show the rms noise amplitude /i„(/) for 


our model (0) of the initial and advanced LIGO inter¬ 
ferometer noise curves, together with the characteristic 
amplitude /ichar(/) for two different examples of BBH co¬ 
alescences, namely a coalescence of total mass 20 Mq at 
a distance of D — 200 Mpc, and a 30 Mq coalescence at 
a redshift of z = 1. (We assume that the cosmological 
parameters are Oq = 1 and ffg = 75 kms“^ Mpc“^.) In 
each case, the sloped portion of the dashed /ichar line is 
the inspiral signal, the flat portion is our crude model of 
the merger, and the separate dotted portion is the ring- 
down. Note that the ringdown and merger overlap in the 
frequency domain since (as we have defined them) they 
are disjoint in the time domain. By contrast, the inspi¬ 
ral and merger are approximately disjoint in both the 
frequency and time domains, since the adiabatic approx¬ 
imation is just breaking down at the end of the inspiral 


Notice that in both cases, 20 Mq and 30 Mq, the 
waves’ characteristic amplitude /ichar(/) is rather larger 
than hn{f) for most of the merger spectrum for the ad¬ 
vanced interferometers, indicating the detectability of the 
merger waveform when optimal filtering can be used. In 
particular, note that the waves should be quite visible to 
the advanced interferometers for the 30 Mq binary even 
though it is at a cosmological distance. Even if such bi¬ 
naries are rare, they are visible to such great distances 
that measurements can be made of a huge volume of the 
Universe, making them a potentially very important and 
interesting source. Cosmological binaries have an en¬ 
hanced SNR in part because the cosmological redshift 
moves their frequency spectrum down closer to LIGO’s 
optimal band; this redshift enhancement effect can some¬ 
times be so strong that for two identical mergers at dif¬ 
ferent distances, the more distant merger can have the 
larger SNR. 

Figure also shows that, of these two example BBH 
coalescences, only the nearby one at a distance of D = 
200 Mpc would be detectable by the initial interferom¬ 
eters. Such coalescences may yield an interesting event 
rate for the initial interferometers—as discussed in the 
Introduction, a coalescence rate of very roughly 3yr“^ 
is expected for BBH coalescences in this mass range 
(M < 30Mq) inside 200 Mpc ||^,|^, although there is 
considerable uncertainty in this estimate. 

A qualitatively different, possibly important type of 
source for the initial LIGO interferometers (and also for 
the advanced interferometers) is the coalescence of black 
hole binaries with masses of order M ^ 100 Mq, as we 
have discussed in the Introduction. (We call such bi¬ 
naries intermediate mass BBHs, as distinct from solar 
mass or supermassive binaries.) In Fig. ^ we show the 
characteristic amplitude /ichar(/) for a hypothetical BBH 
coalescence of total mass 100 Mq at a redshift of z = 0.5, 
corresponding to a luminosity distance of 13 = 2.2 Gpc. 
Note in particular that the initial LIGO interferometer 
noise curve has best sensitivity at ~ 200 Hz just where 
the (redshifted) ringdown frequency is located. We dis¬ 
cuss further in Sec. VI the range and possible detection 
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rates for the initial LIGO interferometers for this type of 
source. 



FIG. 2. A merger of a binary consisting of two 50 Mq black 
holes at a redshift of 2 = 0.5, together with the rms noise 
amplitudes hn{f) for both the initial and advanced interfer¬ 
ometer noise curves for LIGO (see caption of Fig. ^). The 
signal to noise ratios for the inspiral, merger, and ringdown 
stages are about 0, 1.7 and 1.0 respectively for the initial in¬ 
terferometer noise level, and about 11, 52 and 11 respectively 
for advanced interferometers. 


Turn, now, to the detection of supermassive BBH sig¬ 
nals by the possible future space-based detector LISA 
[^,^. LISA can study BBH mergers with far higher ac¬ 
curacy and resolution than the ground based interferom¬ 
eters, because the SNR values are typically much higher 
(> 10^). Our calculation of the inspiral SNR values for 
LISA differs from the other SNR calculations i n the fol¬ 
lowing way. In the frequency integral of Eq. (2.34), if 
one were to integrate over the whole frequency domain 
allowed by the noise curve (as is done for the LIGO initial 
and advanced interferometers), one would in some cases 
obtain the SNR for a measurement which lasts several 
hundred years, which is unrealistic. Thus, it is neces¬ 
sary when calculating inspiral SNR s for LISA to restrict 
the integral over frequency in Eq. (2.34) to the domain 
that corresponds to, say, one year of observation; see Ap¬ 
pendix A. 

Figure^ shows our approximate model [Eqs. ( |4.1| ) and 
(4.4)] of EISA’s projected noise spectrum, together with 
the gravitational-wave amplitude /ichar(/) for the inspi¬ 
ral, merger and ringdown stages of two different BBH 
coalescences: a BBH of total mass IO^Mq at a redshift 
of 2 = 5, and a BBH of total mass 5 x IO'^Mq at a redshift 
oi z = 1. The 10 ®Mq BBH enters the LISA waveband 
a.t f = fs — 10“^ Hz roughly one week before the final 


merger. The SNRs obtained in this case from the inspi¬ 
ral, merger and ringdown signals are approximately 1800, 
4600 and 1700 respectively. The 5 x IO'^Mq BBH enters 
the LISA waveband about twenty years before the final 
merger. The SNR obtained from the last year of the in¬ 
spiral signal, from / ~ 1.6 x 10“^ Hz to / ~ 4 x 10“^ Hz is 
approximately 900, while the merger and ringdown SNRs 
are about 70 and 4 respectively. 



frequency f (Hz) 


FIG. 3. The noise spectrum hn{f) of the possible future 
space-based detector LISA (Laser Interferometer Space An¬ 
tenna), together with the characteristic amplitude /ichar for 
the various stages of two different mual-mass binary black 
hole coalescences (see caption of Fig. nl). The first coalescence 
is that of a binary of total mass 10® Mq at a redshift of 2 = 5. 
The inspiral signal of this binary enters the LISA waveband at 
/ ~ 10““^ Hz about one week before the final merger, and the 
signal-to-noise ratios from the inspiral, merger and ringdown 
are about 1800, 4600 and 1700 respectively. The second coa¬ 
lescence is that of a binary of total mass 5 x 10^ Mq at a red¬ 
shift of 2 = 1, which enters the LISA waveband about twenty 
years before the final merger. In this case a signal-to-noise ra¬ 
tio of approximately 900 would be obtained for the last year 
of inspiral (from / ~ 1.6 x 10““* Hz to / ~ 4 x 10“^ Hz). The 
signal-to-noise ratios from the merger and ringdown in this 
case would be about 70 and 4. 


B. The general signal-to-noise ratio results 

We now turn from these specific examples to the de¬ 
pendence of the SNR values on the mass of and distance 
to the binary in general. In Appendix A we obtain an¬ 
alytic formulae for the SNR values for the three phases 
of BBH coalescences, and for the various interferome¬ 
ters. In this section we plot the results for equal-mass 
BBHs, which are shown in Figs. 11 and The in- 
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spiral and merger curves in these figures (except for the 
LISA insp iral c urves ; see Appendix A) are obtained from 
Eqs. (A4) and (AlO) of Appendix A, while the ring down 
curves are ob tained by numerically integrating Eq. ( 3.24 ) 
in Eq. (pl). 



FIG. 4. The signal to noise ratio (SNR) for equal-mass 
black hole-black hole coalescences detected by LIGO initial 
interferometers, assuming Wiener optimal filtering, as a func¬ 
tion of the redshifted mass (1 -|- z)M of the final black hole, at 
a luminosity distance of 1 Gpc. For fixed redshifted mass, the 
SNR values are inversely proportional to luminosity distance. 
The solid, dotted, and dashed curves are the SNR values from 
the inspiral, merger and ringdown portions of the signal re¬ 
spectively. For non equal-mass binaries, the inspiral SNRs will 
be reduced by the factor ~ ^J4:|l/M, while the merger and 
ringdown SNRs will be reduced by ~ 4p/M; thus the inspiral 
signal will be enhanced relative to the merger and ringdown 
signals. This plot indicates that it is possible that an impor¬ 
tant source for the initial LIGO interferometers may be the 
coalescences of binary black holes where the black holes have 
masses in the range of several hundred solar masses, which 
would be visible out to almost 1 Gpc. For such sources, the 
inspiral portion of the signal would not be detectable, and the 
waves would have to be detected using either the ringdown or 
merger portions. 


The SNR values for the initial LIGO interferometers 
are shown in Fig. |[ This figure shows that it is possible 
that an important source for the initial LIGO interfer¬ 
ometers may be the coalescences of binary black holes 
where the black holes have masses in the range of several 
hundred solar masses, which would be visible out to al¬ 
most 1 Gpc. For such sources, the inspiral portion of the 
signal would not be detectable, and one would need to 
search for the ringdown or merger portions of the signal 
in order to detect the waves. See Sec. VI for further dis¬ 


cussion. The event rate for such high mass BBHs is very 
uncertain; see Ref. for a possible formation scenario. 

Intermediate mass BBHs with /j, <C M {e.g., mi ~ 
10 Mq, m 2 500Mq) are presumably much more com¬ 
mon in the Universe than the intermediate mass BBHs 
with pL ^ M discussed above. The SNRs for such 
mixed binaries will be much lower, however. As seen 
in Appendix A, the merger and ringdown SNR’s scale as 
(/i/M)^, while the inspiral ringdown scales as n/M. The 
difference in scaling occurs because the duration of the 
merger and ringdown signals are approximately indepen¬ 
dent of /i, but the duration of the inspiral signal scales 
as l//r. 

For BBH coalescences of total mass less than or of or¬ 
der 20 Mq, the results of Fig. ^ can be summarized by 
the following approximate for m ulae f or the initia l LIGO 
interferometers [c.f. Eqs. (^), ( AlO ), and ( A16 ) of Ap¬ 
pendix A]: 


/ 5 \ ^ 200Mpc \ ( {l + z)M \ 

UAnspirar ■ V D{z) A 18 Me ) 

(S\ /e^y/2 /200Mpc\ /(l + z)My/" 

■ U.J V D{z) ) 18Me ) 

(S\ / e, n 1/2 p00Mpc\ /(l + z)My/2 

UAingdown^ ■ ^0-032' V D{z) )\ 18 Mg ) ’ 

(5.4) 


where is the fraction of the total mass-energy radiated 
in the merger signal and is the corresponding fraction 
for the rii^down signal. 

Figure 0 also shows that the inspiral signals from low 
mass BBH mergers with M ^ SQMq should be visible 
to ~ 200Mpc (the SNR detection threshold is ~ 5 p0|). 
The ground-based interferometers will, over a period of 
years, gradually be improved from the initial sensitiv¬ 
ity levels to the advanced sensitivity levels |^. Roughly 
half way between the initial and advanced interferome¬ 
ters, the range of the detector system for M < 3QMq 
BBHs will be ^ 1 Gpc. As explained in the Introduction, 
it is thought to be plausible that most BBH progenitor 
systems do not disrupt during the stellar collapses that 
produce the black holes, so that their coalescence rate 
could be about the same as the birth rate for their pro¬ 
genitors, about 1/100,000 years in our Galaxy, or several 
per year within a distance of 200 Mpc . If 

this is the case, BBHs should be detected early in the 
gradual process of interferometer improvement. 

It is also conceivable that BBHs could be seen closer 
than 200 Mpc. High mass stars may have been much 
more numerous among population HI stars due to the 
lower metallicity at high redshift, so there could be a 
large population of BBHs that are remnants of popu¬ 
lation HI stars living in galactic halos. Galactic nuclei 
could also produce a larger population of BBHs than nor¬ 
mal stellar evolution in non-nuclear regions Either 


23 



















of these scenarios could increase event detection rate and 
hasten the first detections of BBHs. 



(l-z)M measured by Advanced LIGO Interferometer (Mg) 


FIG. 5. The SNR values for advanced LIGO interferome¬ 
ters for the inspiral (solid line), merger (dotted line) and ring- 
down (dashed line) phases of equal-mass BBH coalescences; 
see caption of Fig. W For hxed redshifted mass, the SNR 
values are inversely proportional to luminosity distance. For 
values of the redshifted hnal mass lower than ~ 60 Mq the in¬ 
spiral SNR is largest, while for larger BBH systems the merger 
and/or ringdown portions of the signal dominate. 


Figure || shows the SNR values for the advanced LIGO 
interferometers. For BBH coalescences of total mass less 
than or on the order of a few tens of solar masses, the 
results in Fig. can be summarized by the following ap¬ 
proximate formulae for the advanced LIGO interferome¬ 
ters [c.f. Eqs. (|A4D, (AlC) and ( A16| ) of Appendix A]: 



Our estimate (5.5) for the ringdown SNR is roughly con¬ 
sistent with previous estimates by Finn [ p^ . 

From Fig. || it can be seen that for advanced LIGO 
interferometers, equal-mass BBH inspirals will be vis¬ 
ible out to z ^ 1/2 for the entire range of masses 


IOM 0 ^ (1 + z)M < 3OOM0. Thus, there is likely to 
be an interesting event rate. Indeed, the SNRs will be 
high enough even for rather large distances that it should 
be possible to extract each binary’s parameters with rea¬ 
sonable accuracy 0. By contrast, the ringdown SNR is 
fairly small except for the largest mass systems. For very 
massive binaries or binaries that are closer than 1 Gpc, 
advanced interferometers may measure some fairly good 
ringdown SNRs, which would allow fairly good estimates 
of the mass and spin of the final black hole . 

Figure || shows the SNR values obtainable from the 
three phases of BBH coalescences by LISA: the last 
year of inspiral, the merger and the ringdown. We also 
show the SNR value obtainable from one year of inte¬ 
gration of the inspiral signal one hundred years before 
the merger, and a similar curve for one thousand years 
before the merger. This figure shows that LISA will 
be able to perform very high accuracy measurements of 
BBH mergers (SNR values ^ 10^) essentially through¬ 
out the observable Universe {z < 10) in the mass range 
Itf Mq ^ M < 10^ Mq. As discussed in the Intro¬ 
duction, there is likely to be an interesting event detec¬ 
tion rate. The SNR curves in Fig. ^ for measurements 
one hundred and one thousand years before merger show 
that many inspiraling BBHs that are far from merger 
should be detectable by LISA as well. If the merger rate 
of SMBH binaries turns out to be about one per year 
throughout the observable Universe, then at any given 
time one would expect roughly one thousand SMBH bi¬ 
naries to be a thousand years or less away from merger. 
LISA will be able to monitor the inspiral of such binaries 
(if they are of sufficiently low mass) with moderate to 
large SNR | p^ . 

Finally, it should be noted that the relative magnitude 
of the merger and ringdown SNR values is s omewhat un¬ 
certain. As we have discussed in Sec. HI, we have as¬ 
sumed a total radiated energy of O.IM in the merger por¬ 
tion of the signal, and 0.03M in the ringdown portion, a 
ratio of 3 : 1. In reality, it may turn out that in individual 
cases the ratio is as high as 10 or as low as < I. It may 
even turn out to be the case that for most coalescences, 
the ringdown portion of the waveform carries most of the 
radiated energy of the combined merger/ringdown regime 
(depending possibly on the distribution of initial spins). 
Thus, the SNR values shown in Figs. 11 and ^ should 
merely be taken as illustrative. 

We now turn to a discussion of the implications of the 
above signal-to-noise ratios for the required search strate¬ 
gies for the various types of signals, and for the distances 
to which these types of signals can be detected. 
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(l+z)M measured by LISA (Mg) 


FIG. 6. The signal to noise ratio (SNR) for equal-mass 
black hole-black hole coalescences detected by LISA, assuming 
Wiener optimal filtering, as a function of the redshifted mass 
(1 -L z)M of the final black hole, at a luminosity distance of 
1 Gpc. For fixed redshifted mass, the SNR values are inversely 
proportional to luminosity distance. The dotted and dashed 
curves are the SNR values from the merger and ringdown por¬ 
tions of the signal, respectively. The upper solid curve is the 
SNR value that would be obtained from measuring the last 
year of the inspiral signal. For values of the redshifted final 
mass lower than about 10® Mq, the last-year-inspiral SNR is 
largest, while for larger BBH systems the merger and/or ring- 
down portions of the signal dominate. Also shown (lower solid 
curves) are the SNRs that would be obtained from one year 
of integration of the inspiral signal at one hundred and one 
thousand years before the final merger. If the rate of SMBH 
coalescences within 2 < (a few) is roughly one per year, then 
one would expect roughly one thousand SMBH binaries to be 
a thousand years or less away from merger. This plot shows 
that LISA will be able to measure the inspiral of such binaries 
(provided they are of sufficiently low mass) with moderate to 
large SNR M. 


VI. IMPLICATIONS FOR DETECTABILITY OF 
THE GRAVITATIONAL-WAVE SIGNAL 

One of the reasons that coalescences of compact ob¬ 
jects are such good sources for gravitational-wave detec¬ 
tors is that the inspiral portion of the signal is very pre¬ 
dictable, so that optimal filtering may be used in search¬ 
ing for the waves 0. As we have discussed, optimal 
filtering enhances the achievable SNR values by a factor 
of roughly where Afcyc is the number of cycles 


in the waveform in the frequency band of the detectors. 
For neutron star-neutron star (NS-NS) coalescences, Afcyc 
will be on the order of several thousand, while for BBH 
coalescences it will be on the order of several hundred 
[^ . Thus, for NS-NS coalescences, and also for BBH 
coalescences with total mass a few tens of solar masses 
or less, the inspiral portion of the gravitational-wave sig¬ 
nal will be used to detect the entire waveform. In these 
cases, it is not necessary to search for the merger and 
ringdown portions of the waveform, since it will be known 
roughly where in the interferometer data stream they are 
expected to lie. However, there is no guarantee in these 
cases that the merger and ringdown portions will be vis¬ 
ible in the noisy data stream. A key issue is how visible 
these portions are and how much information can be ex¬ 
tracted from them. We discuss this issue in detail in 
paper H [|§. 

For larger mass black hole systems, however, our sig¬ 
nal to noise ratio calculations show that the merger and 
ringdown SNR values can be larger than the inspiral SNR 
values. For equal-mass BBHs, this will occur whenever 
(1 -L z)M > 30Mq for the initial LIGO interferometers, 
and whenever (1 + z)M > 60Mq for the advanced LIGO 
interferometers. Moreover, the inspiral SNR completely 
shuts off for large enough (1 + z)M, as can be seen from 
Figs. ^ and |^. The reason for this shutoff is that at large 
enough values of {l+z)M, the entire inspiral energy spec¬ 
trum lies at frequencies below the “seismic wall” shutoff 
in the interferometer noise curve. The mass cutoff occurs 
at {l + z)M ~ IOOMq for the initial interferometers, and 
at (1 + z)M ~ 400Mq for the advanced interferome¬ 
ters. Admittedly, BBH binaries of total mass ^ 20Mq 
may well be very much more rare than BBH binaries 
of ^ 20Mq; however, they will be visible to such great 
distances that there may be an interesting event rate. 
Moreover, for the initial LIGO interferometers, the mass 
scale ~ 30Mq at which the inspiral SNR becomes much 
smaller than the merger/ringdown SNRs is not terribly 
high. 

In such high mass cases for which the merger and in¬ 
spiral SNRs exceed the inspiral SNR, it will be neces¬ 
sary to perform a search for the merger and/or ringdown 
portions of the signal, independently of any searches for 
inspiral signals, in order that all possible events be de¬ 
tected. If one seeks to detect the waves merely by op¬ 
timal filtering for the inspiral waveform, some fraction 
of the events will be missed which otherwise might have 
been detectable. In fact, it may very well turn out that 
merger signals from BBH coalescences could be the dom¬ 
inant source for the initial LIGO interferometers. 

Therefore a search in real time should be performed 
in the data stream for ringdown and merger signals from 
high mass, BBH mergers. One might imagine that the 
gravitational waves would generally be easier to detect 
by searching for the merger signal, since we have esti¬ 
mated that the SNR values for the merger phase are typ¬ 
ically a factor of a few larger than those for the ringdown 
(c/. Figs. ^ and ||). There are several factors that com- 


25 











plicate this conclusion, however. On the one hand, the 
ringdown’s waveform shape is better understood, which 
makes it easier to produce search templates and hence 
easier to detect the signal. On the other hand, the ratio 
between the merger and ringdo wn S NRs is really quite 
uncertain, as discussed in Sec. V_B, and so it is plau¬ 
sible that the merger SNR will be larger than we have 
indicated relative to the ringdown SNR. In any case, the 
ratio between merger and ringdown SNRs will presum¬ 
ably vary a lot from event to event. Thus, it would 
seem that searches will be necessary for both types of 
signal in the data steam, at least for the mass range in 
which the ringdown SNR is expected to exceed the in¬ 
spiral SNR. [From Sec. ^ we estimate this mass range to 
be {l-\-z)M > 200 Mq for the advanced interferometers, 
and (1 -I- z)M > 60Mq for the initial interferometers.] 
We summarize the discussion of this section by display¬ 
ing the optimum search strategies for various mass ranges 
for the three different interferometers. In each case be¬ 
low, the mass range marked merger refers to matched 
filtering searches for merger signals. If merger templates 
are available, then in the indicated mass ranges merger 
searches will probably be more successful than inspiral or 
ringdown searches; the question mark is a reminder that 
templates for the merger may not be available, and that 
consequently such searches might not be possible. 


ringdown frequency and damping time, can be integrated 
against the (suitably pre-filtered) interferometer output. 
This will be feasible for both the initial and advanced 
LIGO interferometers, and also for LISA. 

The number A/tempiates of required templates can 
be estimated by combining the formalism developed by 
Owen and the results of Echeverria and Finn on the 
expected accuracy of measurement of of the ringdown 
frequency and damping time [p^ , p^ . Using Eqs. (4.15) 
of Ref. |1U| and Eqs. (2.23) and (2.28) of Ref. [^ we 
find that the metric defined by Owen on the space of 
parameters is given by 




(6.4) 
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( 6 . 1 ) 


where Q is the quality factor. The formula ( |6.4[) for the 
Owen metric is valid only in the high Q limit; it has cor¬ 
rections of order IjQ^. Moreover, the formula is also only 
valid when the noise spectrum Sh{f) does not vary sig¬ 
nificantly within the resonance bandwidth A/ ~ /qnr/Q- 
Therefore estimates obtained from Eq. ( |6.4D for the num¬ 
ber of template shapes required for ringdown searches 
will only be accurate to within factors of a few, but this 
is adequate for our purposes. 

Using Eq. (2.16) of Ref. we find that that the num¬ 
ber of required templates is approximately 


A4 


templates 


iQmax(l-MM)-lln 


M„ 


Mrr 


(6.5) 


IMq < M < 200Mq 
200Mq < M < 3000Mo 
SOMq < M < 3000Mq 


IO^Mq < M < IO^Mq 
IO^Mq < M < 10®Mq 
2 X 10®Mq < M < IO'^Mq 


where Qmaxj Mmin and M^ax are the extremal values of 
the quality factor and of the black hole mass that de¬ 
fine the range of signal searches. The quantity MM in 
the formula ( ll.5| ) is the minimal match parameter intro- 
ADVANCEDj^^(,g(j Owen; its meaning is that the event detection 


LIGO 


INTERF 

( 6 . 2 ) 


LISA 

INTERF. 

(6.3) 


rate obtainable by a lattice of templates that have a cer¬ 
tain minimal match MM is smaller than the ideal event 
detection rate (obtainable from an infinitely dense lat¬ 
tice of templates) by the factor (MM)^ |^. We assume 
MM = 0.97 as in Ref. |Q, which corresponds to a 10% 
loss in event rate due to the coarse gridding of the tem¬ 
plate space, and take Qmax = 100 [which by Eq. ( 3.23| ) 
corresponds to 1—a ~ 10“"^]. For the initial and advanced 
LIGO interferometers, the mass range to be searched cor¬ 
responds to roughly Mmin — IMq and M^ax = 5000Mq; 


this yields from Eq. (6.5) 


A. The detectability of high mass black-hole 
coalescences via the ringdown signal 

Consider first the search for ringdown signals. In this 
case, since the shape of the signal is known up to several 
unknown parameters, it will be feasible to implement a 
search based on the method of matched filters, analogous 
to what is being planned for detecting mergers of NS- 
NS and low mass NS-BH and BBH binaries using the 
inspiral waveform p0| , ^ , ^6||^ ] . In other words, a bank 
of templates of the form (3.21), for different values of the 


Ateinplates ^ 4000. 


( 6 . 6 ) 


This is a rather small number of template shapes com¬ 
pared to the number expected to be necessary for inspiral 
searches and therefore the ringdown search should 
be fairly easy to implement. A similarly small number of 
required template shapes (A/tempiates ^ 6000) is obtained 
for LISA assuming M„iin ~ 10^ Mq and M^ax ~ 10®Mq. 

We next discuss the distance to which BBH mergers 
should be dete ctab le via their ringdown signals. As ex¬ 
plained in Sec. IIB, a rough estimate of the appropriate 
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signal-to-noise ratio threshold for detection using one in¬ 
terferometer is @ 


Pthreshold ~ ln [A/templat es T/(eAt)] (6.7) 


where T is the observation time, At is the sampling time 
and e = 10“^ is as defined in Section |ll B| . In fact coin- 
cidencing between the 4 different interferometers in the 
LIGO/VIRGO network will be carried out, in order to 
increase detect ion reliability and combat non-Gaussian 
noise (see Sec. 11 A ). If the noise were exactly Gaussian, 
the appropriate detection criterion would be to demand 
that 


^ ^ P'i — /^threshold’ (^■^) 

3 


fractional loss in event detection rate, if searches for the 
merger signal are not carried out, is the cube of the ratio 
of the SNR values). 

Consider first the ideal situation in which theoretical 
template waveforms are available, so that matched filter¬ 
ing can be used in searches. From Figs, ^and ||it can be 
seen that the merger SNR values are larger than the inspi- 
ral/ringdown values by a factor of up to ~ 4, in the mass 
ranges 30Mq ^ M < 200Mq for initial LIGO interfer¬ 
ometers and IOOMq ^ M < 4OOM0 for advanced LIGO 
interferometers. More precisely, in this mass range, 



( 6 . 10 ) 


where the sum is over the different SNRs obtained in 
each interferometer. In order to combat non-Gaussian 
noise, the detection criterion will be modified to require 
approximately equal SNRs in each interferometer: 

Pj ^ ^threshold/ for all _). (^-l^) 

We have chosen a factor of >72 here to be conservative; 
it corresponds to combining the outputs of just two in¬ 
terferometers (say, the two LIGO 4km interferometers) 
instead of four interferometers. 

Taking T = lO^s and At = 1 ms yields the estimate 
Pthreshoid/v^ « 6.0 for the initial and advanced LIGO in¬ 
terferometers. Therefore, from Fig. we see that the ini¬ 
tial LIGO interferometers should be able to see ringdowns 
from equal-mass BBHs in the mass range IOOMq ^ M < 
700Mq out to ^ 200 Mpc, if the radiation efficiency 
is as large as we have estimated. The advanced LIGO 
interferometers, by contrast, should see ringdowns in the 
mass range 2 OOM 0 < (1 -I- z)M < SOOOMq out to z ~ 1 
(from Fig. For non equal-mass BBHs, these distances 
are reduced by the factor ~ (4/r/M). 

For LISA, the detection threshold is given by Eq. (13 . 
Although LISA does incorporate several partially inde¬ 
pendent interferometers, we have used the noise spectrum 
( [4.4D which is the effective noise spectrum that applies to 
the LISA detector as a whole 0. Thus it is consistent 
to treat LISA as one interferometer. Taking T = 10^ s 
and At = Is, and using the value Afshapes = 6000 
estimated above yields ^threshold ~ 7.5. Hence, from 
Fig. LISA should see ringdowns in the mass range 
IO^Mq < (1 -b z)M < 3 X 10®Mo out to z > 100. 

B. The detectability of high mass black hole 
coalescences via the merger signal 

We next discuss the feasibility of searches for the 
merger signal. As we have explained, this will be most 
necessary when the merger SNR is larger than both the 
inspiral and ringdown SNRs by factors of a few (since the 


The detection threshold for merger searches should be 
approximately the same as that for inspiral and merger 
searches, if the number of template shapes Afshapes is not 
too large (see further discussion below). Therefore, the 
gain in event rate over inspiral/ringdown searches, due to 
searching for merger waves, should vary between 1 and 
about 4^ = 64, depending on the mass of the system, 
if our estimates of and e^. are approximately correct. 
Glearly, the large possible gain in event r ate sh ows the im¬ 
portance of searching for merger waves [1^|. Moreover, 
the ratio of em/^r of energies radiated in the merger and 
ringdown stages, which we have estimated to be roughly 
3, is quite uncertain; we think that in some cases it could 
be rather larger than this (notwithstanding the fact that 
there has been at least one wagered bottl e of wine in the 
community that Cm/cr ^ 1 generically |1^]). In such 
cases the benefit of searching for the merger waves would 
be even greater. 

It is not clear, however, how feasible it will be to pro¬ 
duce a set of numerically generated templates that is 
complete enough to be used to successfully implement 
an optimal filtering search. There may be a very large 
number of distinct waveform shapes, each of which will 
require extensive numerical computations. If both black 
holes are spinning rapidly, the waveforms could depend in 
significant and nontrivial ways on 6 distinct angular pa¬ 
rameters, suggesting that the number of distinct shapes 
could be very large. 

Turn next to the situation where templates are unavail¬ 
able. Gonsider first band-pass filt ering searches. From 
the estimate \/A4.ins = 4 of Sec. HI E| , combined with 
Eq. ( 2.14 ), we obtain the relation 



merger, band—pass search 


1 

4 



merger, optimal 


( 6 . 11 ) 


Therefore, from Eq. ( |6.10D it is apparent that the maxi¬ 
mum achievable band-pass filtering SNR for the merger is 
likely to be essentially no larger than the SNRs obtained 
from the inspiral/ringdown signals. 
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However, as explained in Secs. ID and IIB, noise 


monitoring searches will be more efficient than band¬ 
pass filtering searches, and will be almost as efficient as 
matched filtering searches in most cases, when search¬ 
ing for merger waves. [For the inspiral waves, by con¬ 
trast, noise-monitoring searches perform considerably 
worse than matched filtering searches, since the parame¬ 
ter A/bins is much larger for inspiral waves 1000) than 
it is for merger waves (< 60)]. The event-detection rate 
from noise-monitoring is a factor 


7 ^ = 


f P 




\ /^threshold / 


( 6 . 12 ) 


lower than the event rate from matched filtering. Here 
p* is the noise-monitoring detection threshold, given by 
Eqs. ( 2.30D and (2.31) as a function of the parameters 
e, A^start—times and Adkins; and pthreshoid Is the matched 
filtering threshold, given by Eq. ( ^.15| ) as a function of 
the parameters A'^hapes and A^tart—times- 

Equation ( |6.12| ) is based on the assumption that the 
statistical properties of the noise are Gaussian. In real¬ 
ity the noise in the individual interferometers will have 
significant non-Gaussian components. Coincidencing be¬ 
tween several interferometers is expected to reduce the 
effect of these non-Gaussian components, and the effi¬ 
ciency of this reduction may be higher for matched fil¬ 
tering than for noise-monitoring for the reason discussed 
in Sec. H B . Hence, if coincidencing does not completely 
efface non-Gaussian effects, the loss factor in e vent rate 
TZ may be somewhat higher than the estimate (6.12). 

We now estimate the loss factor in event rate TZ. To 
obtain the most pessimistic estimate possible of the loss 
in event rate, we use the following assumptions: (i) 
The number of template shapes in the matched filter¬ 
ing search is Ayshapes = 1; a more realistic larger number 
would yield a smaller TZ. (ii) The number of frequency 
bins is A/ bins = 60, twice the upper limit estimated in 
Sec. HI E. (iii) The number of starting times in the data 
stream is Aytart-times = 10®, corresponding to a sampling 
time of 0.1s in a data set of one third of a year. Such 
a large sampling time would only be appropriate for the 
largest mass black holes; more realistic sampling times 
will be smaller. Larger values of A4tart-times give smaller 
value s of TZ. (iv) The parameter e in Eqs. (2.30) and 
(|2.15) is e = 10“®. With these assumptions we obtain 


^threshold = 6.8, /9* = 10.3 and the resulting loss in event 
rate factor is 



FIG. 7. This plot shows that factor by which the event 
detection rate is increased when a matched filtering search 
using theoretical template waveforms is carried out, for grav¬ 
itational waves from the merger phase of binary black hole 
coalescences. It is assumed that in the absence of theoretical 
templates, a search is carried out using the noise-monitoring, 
nonlinear filtering method described in the text. Plotted on 
the horizontal axis is the number A/bins = ST A/ of indepen¬ 
dent frequency bins characterizing the space of signals that 
one searches for; T is the maximum expected signal dura¬ 
tion and A/ is the frequency bandwidth. The gain factor 
TZ in event rate, which is plotted on the vertical axis, de¬ 
pends on the number of statistically independent waveform 
shapes in the set of signals one is searching for, which is cur¬ 
rently unknown. This number of waveform shapes can be 
characterized by the effective dimension Nd of the manifold 
of signals, as explained in the text. The upper, solid line 
shows the gain factor in the limit when the number of tem¬ 
plate shapes is small {Md = 0), and is an upper limit on the 
gain factor obtainable from matched filtering. The lower two 
dashed lines show the gain factor when Md = 5 and Md = 10. 
Our best-guess estimate of the number of frequency bins is 
A/bins = 30, corresponding to T = 50M and A/ = I/ IttM) . 
and A/bins is unlikely to be much larger than 10 0 (Se c. ). 
This plot ca n be generated by combining Eqs. (2.3C), (2.31), 
( ^.15| ) and (3.13) of the text, with the assumed parameter 
values ACtart-times = 10® and e = 10“®. 


TZ = 3.5. 

Hence, the availability of templates could increase the 
event detection rate by a factor of at most 4, and more 
typically only 2 over noise-monitoring searches. 


This discussion assumed that the number of waveform 
shapes Ayhapes is small. As this number increases, the 
performance of matched filtering searches gradually wors¬ 
ens until at A4hapes ~ Ayshapes.max, matched filterin g per - 
forms abou t as well as no ise monitoring (see Sec. |H B| ). 
From Eqs. (2.33) and (2.30), the critical value of the num¬ 
ber of shapes is ~ KT^for A/bins = 60, and ~ lO’’ for 
A/iains = 20, assuming A4tart-times = 10®. The number 
of shapes will vary with the signal-to-noise level p. We 
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can define an effective dimension of the manifold of 
signals by the equation 


In [A4hapes(p)] = -^Nd{p) In [l + p^/Ndip)] . (6.13) 


From Eq. ( 2.33| ), we see that Nd{p) is the dimension of 
the equivalent linear space of signals that has the same 
number of distinguishable waveform shapes with SNR 
< p as the true (nonlinear and curved) manifold of merger 
signals. In Fig. ^ we show the gain factor TZ in event rate 
as a function of A^bins, for the values Nd = 0, 5 and 10. 
The true value of Nd is quite uncertain and will vary to 
some extent with p; at high signal-to-noise levels it could 
conceivably be as large as ~ 10. 

By combining the gain factor of 64 discussed above 
with the loss factor TZ, it follows that noise-monitoring 
searches for merger waves should increase the number of 
discovered BBHs by a factor of up to ~ 10 over those 
found from inspiral and ringdown searches. 

For LISA, the expected SNRs are so high that the avail¬ 
ability of merger templates will likely have no impact on 
event detection rates. However, for LISA, as well as for 
ground-based interferometers, accurate merger templates 
will be vitally important for interpreting the detected 
waveforms, for extracting all the available information, 
and for testing general relativity M] . 


VII. CONCLUSIONS 

It seems quite likely that the gravitational waves from 
merging BBH systems will be detected by the ground- 
based interferometers that are now under construction. 
The initial LIGO interferometers will be able to detect 
low mass (^ 30Mq) coalescences of equal-mass BBHs out 
to about 200 Mpc via their inspiral waves, and higher 
mass (IOOMq ^ M < 700Mq) systems out to about 
200 Mpc via their ringdown waves. The advanced LIGO 
interferometers will be able to detect equal-mass BBH 
coalescences in the mass range IOMq < M < 3 OOM0 
to z ~ 1/2 via their inspiral waves, and higher mass 
( 2 OOM 0 % M 3 OOOM 0 ) systems to z ~ 1 via their 
ringdown waves. For non-equal mass BBHs, these dis¬ 
tances will be reduced by a factor of ~ for in¬ 

spiral signals and ~ 4p/M for ringdown signals. 

Searches for massive BBHs {M > 5 OM 0 for 
LIGO/VIRGO) based on merger waves could increase 
the range of the interferometers by a additional factor of 
~ 2, without requiring detailed knowledge of the wave¬ 
form shapes. It seems likely that BBH coalescences will 
be detected early in the gradual process of improvement 
from the first interferometers to the advanced interferom¬ 
eters, rather than later, and there is a strong possibility 
that they will be the first sources of gravitational radia¬ 
tion to be detected. 

Theoretical template waveforms obtained from numer¬ 
ical relativity supercomputer simulations will be crucial 


for analyzing the measured merger waves. A match of the 
detected waveform with a predicted waveform would be 
a triumph for the theory of general relativity and an ab¬ 
solutely unambiguous signature of the existence of black 
holes. A complete set of such theoretical templates would 
also aid the search for BBHs, but not by a large amount. 

The space-based interferometer LISA will be an ex¬ 
tremely high precision instrument for studying the coales¬ 
cences of supermassive BBHs. Coalescences with masses 
in the range 1O®M0 < (1 -f z)M < IO^Mq should be 
detectable out to z ~ 10 with very large SNRs (> 10^), 
via their merger and ringdown waves. Additionally, sys¬ 
tems in the mass range 1O"^M0 < (1 -|-z)M 3 x 1O’’M0 

should be detected to similar distances and with SNRs 
^ 10^ via their inspiral waves. 
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APPENDIX A: SIGNAL TO NOISE RATIO FORMULAE 


In this appendix, we give the details of our calculations of the signal-to-noise ratios. Th e res ults of this appendix 
were used to make Figs. || andH in the text, and also the approximate SNR formulae (^) and (^) given in Sec. 0. 
The following three subsections give the calculations for the inspiral waves, the merger waves and the ringdown waves 


respectively. 

Note that throughout this appendix we use 


‘MMq” (Mega solar-mass) as shorthand for 10® M©. 


1. Inspiral 


To calculate the angle-a vera ged SNR squared for the inspiral, we insert the inspi ral energy spectrum ( 3.2C ) and 
our parameterized model (0) of an interferometer’s noise spectrum into Eq. (|2.34|) , and integrate from / = /s to 
/ = /merge- The result is 


(P^) = 


T,{M,z,D)v^ 

T,{M,z,D) 

T,{M,z,D) 

0 


9^1/3 _ f - |a-i/3„io/3 _ ^3 

9a3/3-8(^)'/"-a3(^) 
„l/3(^^)-„3(A)®/" 


8/3' 


O^fm ^ /merge/(l T ^') 

/m/o; ^ /merge/(l -f 2:) < 0^/77 

Is — /merge/(l “t“ 2) < /m/ 
/merge/(1 + z) K fs, 


(Al) 


where 


and 


(1 -f z)af, 

V = 


T,{M,z,D) = 


/merge 

[{1 + z)M]^/^[Afi/M] 
SOn-^/^Dizykl^f, 


1/3 


(A2) 


(A3) 


Here D(z) is t he luminosity distance to the sour ce, /s, a, fm and hm are parameters characterizing the detector noise 
spectrum ( 14T|) , and /merge is given by Eq. (Q. 

Inserting the values of the noise spectrum parameters from Eq. (4.2) for initial LIGO interferometers, we obtain 
the following numerical values for the SNR in the equal-mass case /i = A//4: 


= 


initial 


2.8 

4.7 

2.7 
0 


( 200 Mp( 

= ( 

{1+z)M\ 

V 

) 1 

18 Mq j 

( 200 Mp{ 

= ( 

(1+z)M\ 

V 

) \ 

18 Mq j 

( 200 Mp( 


il+z)M\ 

V 

) 1 

36 Mq j 


5/6 

- 1/2 


-I 1/2 


1 /n'i/2 

1 _ 0 71 (" (1+^)^^ 

1 18M0 J 

. . . / 3/3 


-I 1/2 


{l + z)M < 18 Mq 

18 Mq < {l + z)M < 36 Mq 

36 Mq<{1 + z)M < 102 Mq 
102 Mq < {l + z)M. 


(A4) 


For the noise curve parameters (4.3) appropriate for advanced LIGO interferometers we obtain 

31 1/2 


= 


advanced 


/ 

27 

( lGpc\ / (1-I-z)M\ 
^ D{z) J 37 Mq j 

43 

( lGpc\ / (1-I-z)M\ 
^ D{z) ) 37 Mq ) 

31 

( lGpc\ ( (1-I-z)M\ 
i, D{z) J {95 Mq j 

0 



5/6 

- 1/2 


1 - 0.16 

< (1+z)M\ 

< 37Mq j 

1 - 0.65 

f (1+z)M\ 
37 Mq j 


10/3 
1/3 
\ 8/3 


1/2 


-I 1/2 


{l + z)M < 81 Mq 

81 Mq < (l-hz)M< 95 Mq 

95 Mq < (1 -h z)M < 410 Mq 
410 Mq < {l + z)M. 


(A5) 


As explained in Sec. V A, the calculation of the inspiral SNR for LISA differs from the other SNR calculations in 
the following way. If one were to integrate over the whole frequency domain in the interferometer waveband up to 
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/ = /merge (as is done for the initial and advanced interferometers in LIGO), in some cases one would obtain the 
SNR for a measurement of several hundr ed ye ars duration, which is obviously irrelevant. Thus, it is necessary to 
restrict the integral over frequency in Eq. (2.34) to the domain that corresponds to, say, one year of observation when 
calculating inspiral LISA SNRs. Using the Newtonian relationship for the rate of frequency sweep, we obtain for the 
frequency at time T before merger in the equal-mass case 


/insp(T) — 




1 -3/8 


(A6) 


Binaries of redshifted total mass {l + z)M larger than about 5 x IO^Mq enter the LISA waveband at / = /^ = 10“^ Hz 
less than one year before merger, while binaries of s maller red shifte d mass spend more than one year in the LISA 
waveband. To calculate the SNR, we insert Eq. ( 3.20|) into Eq. ( 2.34 ) and integrate numerically from the larger of fs 
and /insp(lyr) to /merge- The resulting SNR values are shown in Fig. We also show in Fig. ^ the SNR obtained 
from one year of observation one hundred years before the final merger, obtained by integrating from /insp(lOOyr) to 
/insp(99yr), as well as a similar curve for one thousand years prior to merger. 

Equation ( |Al| ) applies to LISA only for (1 -|- z)M ^ 5 x IO^Mq. By combining Eqs. (|Al|) and ( f4.4[ ) for (1 -|- z)M > 
5 X IO^Mq together with an approximate fit to Fig. for (1 -1- z)M < 10®Mq we obtain for the SNR from the last 
year of inspiral in the equal-mass case 


^ ) LISA 


1 3 X 10^ ( IHEe') / (l+z)M \ 
i.O X iU j (^o.SMMq j 
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O.U X iU ^ J 6MM0 J 


I - 0 
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0.5 MM© J 


1/3 


1 — 0 006 ('l 


1/2 

ni/2 


IOOM 0 < (I -b z)M < 0.5 MM© 
0.5 MMq < (1 -b z)M < 6.0 MM© 

6.0 MM© < (1 -b z)M < 41MM© 
41MM© < (1 -b z)M. 


(A7) 


2. Merger 


To calculate the merger SNR we use the energy spectrum ( |3.14[) and follow the same procedure as above. The 
result is 


ip^) = 


K'^-l 


Tm{em,M,z,D)v^ 

Tm{im,M,z,D) 31ni;- 


^Tm{(lm,M,Z,D) 
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(*) +31"(“f) 


[a^v - 1)] 


Tm{^m,M,Z,D) 
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£/s 


/merge/(l ■2) ^ ^fm 

fm!^ ^ /merge/(l H“ 2 ) C^fm ^ /qiir/(l “1“ 
/merge/(l ^ fmjOifm /qnr/(1 H“ 

fmjOt /merge/(l z) /qnr/(1 H“ ■2^) 

/merge/(l ^ fs ^ fmj^ ^ /qnr/(l 

fs ^ /merge/(l ~\~ z') fm!^ ^ /qnr/(l ~\~ zfj 

/merge/(l 2^) ^ /s fm!^ /qnr/(l -^) 

fs ^ /merge/(l H“ 2 ) <C /qnr/(l H“ 

/merge/(l ^ f s ^ /qnr/(l H“ ^ fmj^ 

/qnr/(l + z) < fs- 


(A8) 


Here v is given by Eq. ([A^), Cm is the fraction of total mass energy radiated during the merger (which we have also 
denoted by Cmerger in the body of the paper), k = /qnr//merge, and 


Truism, M,Z,D) = 


2e^M(l + z)"[4^/M]" 
157r2D(z)^/l^/merge(^^ - 1) 


(A9) 
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Lines marked with the superscript “I” turn out to hold for the initial LIGO interferometer parameters; those with 
“A” hold for advanced LIGO interferometer parameters; and t hos e with “L” hold for LISA. 

U sing the num eric al values of the noise curve parameters (4.2) for initial LIGO interferometers, and Eqs. ([l.4|), 
(3.5), (A2), and (^) we find for the initial LIGO interferometers in the equal-mass case 
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Similarly using Eq. ( |4.3| ) we find for advanced LIGO interferometers 
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620 Mq < (1 -f z)M < 2600 Mq 
2600 Mq < {l + z)M. 


(All) 


Finally, using the parameters (^) appropriate for LISA, we obtain 
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(A12) 


3. Ringdown 


The ringdown S NRs are calculated a little differently from the inspiral and merger SNRs. First, we use the e ffective 
energy spectrum (3.24) which yields an estimate of the true SNR obtainable from the model waveform (3.21) that is 
accurate to within a few tens of percent (see Appendix B). Second, the integral over frequency in the SNR formula 
( ^.34 ) with the noise spectrum ( ^ ) and the energy spectrum (3.24) cannot easily be evaluated analytically. Henc^ 
we calculated this integral numerically to produce the plots of ringdown SNR versus BBH mass shown in Figs. 
and ||. 

In the remainder of this appendix we derive approximate formulae for the ringdown SNR as a function of mass, 
by approximating the ringdown energy spectru m as a delta function at the ringdown frequency [cf. Eq. (3.25)]. This 
approximation yields (see Appendix B and Ref. |103[) 


(P") = 


(1 + z)3mM"Q[4/i/M]3 

20Tr'^D{z)‘^fqnrSh[fqnY/il + z)] ' 


(A13) 


Using Eq. ( 3.23 ) and the relation ( 3.28| ) between the dimensionless coefficient A and the radiated energy we can 
rewrite formula (A13) as 
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8 1 (1 -f z)M 

'{l + z)M' 
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5 F(a)2 "" ^,.[/qnr/(l + ^)] 

[ D{z) \ 


_M_ 


(A14) 


where is the fraction of the total mass energy radiated in the ringdown, and 


^ ' 100^ ' 


(A15) 


An equivalent formula was previously obtained by Finn |102]. 

We find the following numerical result when we insert our assumed values = 0.03 and a = 0.98 for the ringdown 
signal together with the parameters for the initial LIGO interferometer noise curve in the equal-mass case: 
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The corresponding formulae for advanced LIGO interferometers are 


= 


advanced 
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Finally, the corresponding formulae for LISA are 
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(A17) 


(A18) 


By comparing Eqs. (A16) - (A18) with Figs. it can be seen that the delta-function energy spectrum approximation 

is fairly good except for M > 3 OOOM 0 for advanced LIGO interferometers and M > 3 x 10®M 0 for LISA. The 
approximation fails to capture the high mass tails of the SNR curves. 
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APPENDIX B: ENERGY SPECTRUM FOR 
RINGDOWN WAVES 

There is a subtlety in calculating the SNR for the ring- 
down waves, related to the fact that the SNR squared 
does not accumulate locally in the time domain. In or¬ 
der to explain this subtlety, let us focus not on the angle- 
averaged SNR squared which was our main concern in the 
body of the paper, but rather on the SNR squared ob¬ 
tained in one interferometer from a specific source with 
specific relative angular orientations. In this case the 
waveform h(t) seen in the interferometer, for t > 0, is of 
the form 


h{t) = ho cos(27r/q„rt -|- V'o)e 


(Bl) 


for some constants ho and ipo, while h{t) is the (unknown) 
merger waveform for t < 0. 

Let us also focus first on the simple, idealized case 
of white noise, Sh{f) = Sh = const. Then, the SNR 
squared ( 13 ) accumulates locally in time: 


p^ = — dth(tf. 
J — oo 


(B2) 


Hence, for white noise, the SNR squared from the ring- 
down is clearly unambiguously given by 


Z _ 

rringdown ~q 


2Sh 


nOO 

/ hi cos^(27r/qnrt + 
n Jo 

cos{2'ijjo) — Qsin{2'ijjo) 


1-f 


1 + Q2 


§^[1 + 0 ( 1 / Q )] 


(B3) 


where Q = TrfqnrT. Now consider the case when the 
noise is not exactly white. Naively, we expect that in 
the Fourier domain the energy spectrum of the ringdown 
signal will be a resonance curve that peaks at / = /q„r 
with width ~ fqnr/Q- Thus, for large Q we would expect 
that most of the SNR squared will be accumulated near 
/ = /qnr, unless the noise spectrum varies very strongly 
with frequency. Moreover, if the noise spectrum Sh{f) 
does not vary much over the bandwidth ~ fqnr/Q of the 
resonance peak, then we would expect the formula (p33| ) 
to be valid to a good approximation, with Sh replaced 
by Sh{fqni)- We show below that this is indeed the case: 
under such circumstances, Eq. (|B3| ) is fairly accurate, and 
the resulting approximate ringdown SNR is embod ied in 
our approxi mate del ta-fu nction energy sp ectru m ( 3.25| ) 
and in Eqs. ( Alg ) - ( A17 ) of Appendix A 1051. 

In many cases of interest, it will indeed be true that 
most of the SNR squared for ringdown waves will be accu¬ 
mulated in the vicinity of the resonance peak, so that the 
SNR will approximately be given by Eq. (B3). However, 
this will not always be the case. For instance, suppose 
that we were lucky enough that two 10®M© black holes 
were to merge at the center of our own galaxy. Would 


such an event be detectable by advanced LIGO interfer¬ 
ometers? Clearly, most of the power in the ringdown 
waves in this case would be far below the LIGO /VIRGO 
waveband. However, given that the merger is only at 
^ 10 kpc, one might hope to be able to detect the tail of 
the ringdown waves that extends upwards in frequency 
into the LIGO/VIRGO waveband. Or, consider the de¬ 
tectability of a ringdown of a nearby IO^Mq black hole by 
LISA. In this case most of the ringdown power is concen¬ 
trated at frequencies above LISA’s optimum waveband, 
and the detectability of the signal is determined by the 
amount of power in the low frequency tail of the ring- 
down. In such cases, it is clearly necessary to understand 
the power carried in the ringdown waves at frequencies 
far from the resonant frequency. 

Normally, such an understanding is obtained simply by 
taking a Fourier transform of the waveform h{t). In the 
case of ringdown waves from BBH mergers, however, the 
waveform for t < 0 is the unknown, merger waveform. 
In order to obtain the SNR squared from the ringdown 
signal alone, one might guess that the appropriate thing 
to do is to take h{t) = 0 for t < 0, and insert this together 
with Eq. (Bl) into the standard formula (2.7) for the 
signal to noise squared. However, the resulting energy 
spectrum has unrealistic high frequency behavior due to 
a discontinuity in h{t) at t = 0 (or a discontinuity in 
h'(t) at t = 0 in the case ipo = ■^/2), and the resulting 
SNRs can in some cases differ from the correct values 
(see below) by factors > 10. Other choices for h{t) for 
t < 0 [for instance h{t) = /i(0)] get around this problem 
but instead have unrealistic low frequency behavior. In 
any case, it is clear that these choices are somewhat ad 
hoc and that there should be some more fundamental, 
unique way to calculate the SNR. 

We now explain the correct method to calculate the 
SNR. The question that effectively is being asked is: 
What is the probability that there is a ringdown wave¬ 
form present in the data stream, starting at (say) t = 0? 
This probability is to be calculated given only the data 
from t > 0, without using the measured data from t < 0 
which is contaminated by the unknown merger waveform. 
To do this one must effectively integrate over all possible 
realizations of the noise for t < 0. The necessity for such 
an integration is illustrated by the following simple anal¬ 
ogy. Suppose that one is measuring two real variables, 
h+ (“waveform for positive t”) and (“waveform for 
negative t”), and that the measured values of these vari¬ 
ables are h+ and h-. Suppose that because of the noise 
in the measurement process, the probability distribution 
for the true values of these parameters given their mea¬ 
sured values is 


P^h+,h_) = ^ exp I - ^ [{h+ - h+f 

+ {h_ - h_f + 2e{h+ - h+){h_ - /i-)] |. (B4) 
Thus, h+ and /i_ are Gaussian distributed about their 
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means /i+ and /i_, and they are correlated. If we assume 
that h- = 0 [analogous to assuming h{t) = 0 for t < 0], 
we obtain for the probability distribution for 


p{h+\h_ = 0) = 


1 


\f^o 




(B5) 


where h'_^_ = /i+ — eh-. By contrast, if we instead cal¬ 
culate the probability distribution for /i_|_ alone by inte¬ 
grating over h- we find 


p{h+) = 


1 




^-(h+-h+f/{2a-f 


(B6) 


where a* = cr/Vl — It is clear in this sim ple e xample 
that one should use the reduced distribution (^) rather 
than the distribution (^). Note also that the widths of 
the probability distributions a nd (p^ are different, 
and that the correct distribution ((Bq ) could not have 
been obtained from the joint distribution (B4) for any 
assumed choice oi h-. 

Turn now to the analogous situation for random pro¬ 
cesses. If n{t) is the interferometer noise, let Cn{^t) = 
{n{t)n{t + At)) denote the autocorrelation function. De¬ 
fine the inner product 


{hi \ / 12 ) = f 
Jo 


dt / dt' K{t,t')hi{t)h2{t') (B7) 

0 Jo 


on the space of functions h{t) for t > 0, where the kernel 
is determined from 


nOO 

/ df'Kit, - t') = S{t - t') 

Jo 


(B8) 


for t,t' > 0. The quantity is analogous to the 

modified wid th a* in Eq. (|B6| ) above. Using the inner 
product (^), the usual theory of optimal filtering [|^j^ 
can be applied to random processes on the half line t > 0. 
Thus, if s( t) is the interferometer output and h{t) is the 
waveform (BI), the detection statistic is U = {s\h), and 
the SNR squared for the measurement is 




E[Yf 


E[Y^] - £;[y]2 

= {h\h) 


/ dt dt'K{t,t')h{t)h{t'), 
'0 Jo 


(B9) 


where means expectation value. If we define 

pOO pOO 

G{f,f)= dt dt'e^^^f*e-^^^f'*'K{t,t') (BIO) 
^0 Jo 


and 


Hf) = f 

Jo 


h{t) dt, 

the SNR squared can be rewritten as 


(Bll) 


P^ = 


df 


dfHf)*G{f,f)h{f). (B12) 


Note that the Fourier transform G{f, f) of K{t, t') is not 
proportional to 5{f — f')/Sh{f) but instead is in general 
non-diagonal in frequ ency. 

The formula (B12) resolves the ambiguities discussed 
above in the method of calculating the ringdown SNR; 
the result does not require a choice of the wavefo rm h{t) 
for t < 0. Unfortunately, the final answer ( BI2 ) is com- 
plicate d in the sense that it cannot be expressed in the 
form (2.34) for any effective energy spectrum dE/df. 
This is somewhat inconvenient for the purposes of this 
paper: the wave’s energy spectrum is a useful and key 
tool for visualizing and understanding the SNRs. Clearly, 
an approximate, effective energy spectrum (to the extent 
that one exists) would be very useful. We now turn to a 
derivation of such an approx imat e, effective energy spec¬ 
trum, namely the spectrum ( 3.24 ) which is used through¬ 
out the body of this paper. 

We start our derivation by describing an alternative 
meth od o f calc ula ting the exact ringdown SNR given by 
Eqs. (M) and (^^). It is straightforward to show that the 
quantity ( |B9D can be obtained by (i) choosing any wave¬ 
form h{t) for t < 0, (ii) calculating the SNR from the 
usual formula (0), and (iii) minimizing over all choices 
of the function h{t) on the negative real axis. We have 
experimented with several choices of h{—t) for t > 0, 
namely h{—t) = 0, h{—t) = h{0), h{—t) = h{t). We 
found that the SNR obtained by minimizing over these 
choices is always Tor the entire black hole mass ranges 
discussed in Sec. 0) within a few tens of percent of the 
SNR obtained from the following prescription: (i) As¬ 
sume that h(t) for negative t is identical to the waveform 
for positive t except for the sign of t/r; ie., that 


h{t) = ho cos(27r/qnrt + 4’o)e 


(B13) 


for positive and negative t. ( ii) C alculate the total SNR 
using the standard formula (2.7). (iii) Divide by a cor¬ 
rection factor of in amplitude to compensate for the 
doub ling up. This prescription gives the correct, exact re¬ 
sult (M) for white noise. For more realistic noise curves, 
the errors of a few tens of percent resulting from this pre¬ 
scription are unimportant compared to the uncertainty 
in the overall amplitude A of the ringdown signal. More¬ 
over, the resulting SNR values multiplied by are an 
upper bound for the true SNR (since if our ad hoc choice 
of h{t) for t < 0 happened to be exactly right, then the 
prescription would underestimate the SNR by -\/2). 

We now explain how to obtain the energy spectrum 
( |3.24 ) from the above approximate prescription. From 
Eqs. (2.35) and (3.21) it can be seen that the waveform 
as seen in one int erferometer, before angle averaging, is 
given by Eq. (Bl) with 


AM 


(B14) 
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Here the angle s 0, tp, ip, i and /3 have the meanings ex¬ 
plained in Sec. |llQ. L et us now insert the waveform (Bl) 
into the formula (B9) for the exact SNR, and t hen aver¬ 
age ov er the angles 0, ip, t and /3 using Eqs. ( ^.3^ ) and 
( 3.22 ). This yields for the angle-averaged, exact SNR 
squared 


{Pexact[h(t)] ) — [Pexact[h+fi{t)] 


where Pexact[h{t)] denotes the exact SNR functional ( p^ ) 
and 


■ Pexact [^X ,0 (^)] ] ; 

(B15) 


h+,o(t) = cos(27r/qnrt) e 

hx,o{t) = sin(27r/qnrt) 


(B16) 


for t > 0. Now, for e ach of the two terms on the right 
hand side of Eq. (B15), we make the app roxim ation dis¬ 
cussed above consisting of using Eqs. (^) and ( B13 ) and 
dividing by 2. This yields 


{PexactiHt)]^) ~ / df 


|^+,o(/)P + |^+.o(/)P 


ShU) 


(B17) 


where it is understood that /i+,o and /ix,o have been ex¬ 
tended to negative t in the manner of Eq. ( |B13 ). Fi¬ 
nally, evaluating the Fourier transforms yields an angle 
averaged SNR squared of th e form (2.34), with energy 
spectrum given by Eq. (3.24|). 
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yields the formula (A13) for the ringdown SNR, corre- 
spondi ng to the approximate delta function energy spec - 
trum (3.25), proceed as follows. Starting from Eq. ( B14 ), 
avera ge ov er the angles 9,ip,ip,L and /3 using Eqs. (2.3S) 
and (3.22). This yields 


{ho{e,ip,iij,L,pf)= ^ 
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Substitute this into Eq. (B3), make the replacement Sh 
Shifqni), and make the appropriate substitutions to a c- 


count for co smol ogical effects as explained in Sec. HG 
obtain Eq. (A13). 


to 


vWs, as discussed after Eq. ( |2.34| ) above. 
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